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Abstract 

The  linearized  expansion  of  an  Infinitely  long  cylin- 
drical piston  Into  a  conducting  fluid  Is  considered.   The 
surface  of  the  piston  Is  assumed  to  be  a  perfect  conductor. 
The  fluid  Is  assumed  to  be  a  perfect  compressible  fluid 
(free  from  viscosity  and  heat  conduction)  and  a  perfect 
electrical  conductor,  permeated  by  a  uniform  magnetic  field. 

Conical  solutions  of  the  linearized  Lundqulst  equations 
are  obtained  by  an  extension  of  the  conical  plane  wave  method 
first  Introduced  by  Gardner.   In  general,  the  solution 
Involves  both  real  and  complex  characteristics.   In  the  course 
of  the  analysis,  the  correct  way  to  handle  such  coupled  prob- 
lems is  elucidated. 

We  consider  two  sorts  of  problems:   slow  expansions  of 
pistons  with  roughly  circular  cross-section,  and  finite-speed 
expansions  of  thin  pistons.   The  linearized  problem  appears 
to  be  a  faithful  model  of  the  nonlinear  problem  only  in  the 
latter  case,  i.e.,  only  for  pistons  whose  cross-section  is 
nearly  parallel  to  the  undisturbed  magnetic  field.   In  other 
cases,  the  presence  of  a  magnetic  boundary  layer  is  suggested. 
The  linearized  solution  for  a  thin  piston  with  finite  expan- 
sion speed  is  obtained,  and  appears  to  be  physically  relevant. 
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Introduction 

In  this  thesis  we  consider  the  linearized  expansion  of 
an  Infinitely  long  cylindrical  piston  Into  a  conducting 
fluid.   We  assume  that  the  fluid  Is  a  perfect  compressible 
fluid  (no  viscosity  or  heat  conduction)  and  a  perfect 
electrical  conductor,  permeated  by  a  constant  magnetic  field. 
We  obtain  conical  solutions  of  the  linearized  Lundqulst 
equations . 

The  problem  was  suggested  by  two  papers  of  Bernstein, 
Kulsrud,  et  al.  [2,  7].   They  obtained  numerical  solutions  to 
the  three-dimensional  problem  of  an  expanding  sphere.   It  was 
felt  that  their  results  were  inconclusive,  and  that  an  exam- 
ination of  a  two-dimensional  problem  would  shed  light  on  the 
nature  of  the  problem  under  consideration.   We  obtain  conical 
solutions  by  an  extension  of  Gardner's  plane  wave  method  [5]. 

The  problem  is  particularly  interesting  since  there  are 
both  real  and  complex  characteristics.   Separate  representa- 
tions for  the  portions  of  the  solution  which  depend  on  the 
real  and  the  complex  characteristics  are  obtained.   We  first 
consider  slowly  expanding  bodies.   By  solving  a  free-boundary 
problem  we  obtain  a  solution  which  depends  only  upon  the 
complex  characteristics;  the  resulting  piston  cross-section 
is  approximately  circular.   This  solution  suggests  the  presence 
of  a  magnetic  boundary  layer  near  the  piston,  since  the 
magnetic  field  perturbations  near  the  piston  are  of  the  same 
order  as  the  underlying  field.   For  general  reasons  no 
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sol-utlons  Involving  only  the  real  characteristics  can  be 
constr-ucted. 

For  slow  expansions  involving  pistons  of  arbitrary 
cross-section,  both  the  real  and  complex  characteristics 
must  be  employed.   Using  a  non-conf ormal  mapping  of  the 
(^,  1^) -plane  [5]  into  a  complex  z-plane,  we  can  reduce  the 
problem  to  the  determination  of  an  analytic  function  from 
its  boundary  values.   We  present  a  constructive  method  for 
the  solution  of  the  boundary  value  problem. 

Finally,  we  consider  the  finite-speed  expansion  of  a 
thin  piston  which  is  nearly  parallel  to  the  undisturbed 
magnetic  field.   In  this  case  we  reduce  the  problem  to  the 
same  sort  of  boundary  value  problem  mentioned  above.   How- 
ever, in  this  case,  the  boundary  data  is  of  lower  order 
than  the  unperturbed  flow,  so  that  the  resulting  solution 
appears  to  be  physically  relevant.   For  this  problem,  the 
linearized  problem  is  probably  a  faithful  model  of  the 
nonlinear  problem.   Treating  the  piston  as  a  slit,  we  can 
write  down  the  solution  explicitly  in  terms  of  integrals  of 
singular  integrals. 
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1 .   General  Formulation  of  the  Problem 

In  this  first  section  we  shall  discuss  some  of  the 
general  features  of  the  problem  under  consideration,  and 
put  on  record  some  general  properties  of  magnet ohydrodynamic 
flows  which  we  shall  need  later.   For  the  most  part,  this 
section  is  modelled  after  Gardner's  report  [5]>    but  a  good 
many  of  the  missing  details  have  been  supplied. 

The  first  subsection  discusses  the  linearization  of  the 
Lundqulst  equations,  and  the  reduction  to  dependence  on  two 
space  variables.   We  restrict  our  discussion  to  slow  expan- 
sion velocities  of  the  piston  so  that  linearization  is 
permissible.   It  does  not  appear  that  the  linearized  solution 
is  always  a  faithful  model  of  the  nonlinear  solution.   This 
matter  will  be  discussed  in  more  detail  later.   The  Alfve'n 
wave  (given  by  the  components  of  u  and  B  along  the  axis  of 
the  piston)  is  decoupled  from  the  remaining  linearized 
Lundqulst  equations  to  obtain  the  desired  two-dimensional 
formulation . 

Subsection  two  discusses  plane  waves  and  the  construction 
of  the  characteristic  conoid.   We  shall  form  our  solution  by 
superposing  conical  plane  waves.   The  characteristic  conoid 
will  play  a  central  role  in  our  subsequent  discussions, 
since  it  gives  the  range  of  influence  of  a  disturbance  at 
(x,y,t)  =  (0,0,0)  for  the  linearized  Lundqulst  equations. 
Parametric  equations  for  the  characteristic  conoid  are  obtained 


3 


The  third  subsection  discusses  the  representation  of 
a  conical  flow  by  plane  waves.   We  rely  heavily  upon  Gardner's 
work  [5]  In  this  development.   We  obtain  a  formal  representa- 
tion of  a  conical  solution  of  the  linearized  Lundqulst 
equations.   This  representation  forms  the  basis  of  our 
subsequent  solutions  to  the  piston  problem. 

Subsection  four  Is  devoted  to  a  statement  of  the  well- 
known  jump  relations  for  the  Lundqulst  equations  [1].  We 
shall  have  occasion  to  employ  these  relations  In  the  formula- 
tion of  appropriate  boundary  conditions  for  the  problems 
which  we  shall  consider.   They  also  play  a  central  role  In 
the  proof  of  our  uniqueness  theorem. 

The  fifth  subsection  discusses  the  boundary  conditions 
which  we  shall  Impose  at  the  surface  of  the  piston.   Firstly, 
we  assume  that  the  axis  of  the  piston  Is  perpendicular  to  the 
constant  magnetic  field  of  the  undisturbed  flow.   Secondly, 
we  assume  that  the  tangent  to  the  cross-section  of  the 
piston  (at  fixed  time)  Is  perpendicular  to  the  constant 
undisturbed  magnetic  field  at  only  a  finite  number  of  Isolated 
points.   Finally,  we  assume  that  E+uxB=Oon  the  surface 
of  the  piston.   The  first  two  assumptions  make  the  formula- 
tion  of  boundary  conditions  on  the  surface  of  the  piston 
particularly  simple,  and  are  suggested  by  Gardner's  formula- 
tion of  the  half -plane  diffraction  problem  [5].   The  last 
requirement  is  equivalent  to  the  assumption  that  the  surface 
of  the  piston  is  a  perfect  conductor. 


-  4 


The  final  subsection  presents  the  uniqueness  theorem 
for  the  piston  problem.   It  Is  modelled  after  Gardner's 
theorem  [5],  and  employs  the  so-called  abc  method.   Armed 
with  this  uniqueness  theorem,  we  can  be  assured  that  we  are 
proceeding  correctly  so  long  as  we  can  exhibit  a  solution 
satisfying  whatever  (appropriate)  conditions  we  may  prescribe 


I'l .   Equations  of  Motion 

We  assume  that  the  medium  In  which  the  piston  expands 
Is  a  perfect  compressible  fluid  (we  neglect  heat  conduction 
and  viscosity)  and  a  perfect  electrical  conductor.   In  the 

undisturbed  state  the  fluid  Is  at  rest,  and  Is  permeated  by 

_^ 
a  constant  magnetic  field  B  .   We  shall  take  the  axis  of 

*=  o 

the  cylindrical  piston  to  be  perpendicular  to  the  direction 
determined  by  the  undisturbed  magnetic  field,  B  .   The  state 
of  motion  of  such  a  fluid  Is  specified  by  the  following 
quantities : 

p   the  scalar  pressure 

p   the  density 

u   the  fluid  velocity 

B  the  magnetic  field 

E  the  electric  field 

J   the  electric  current  density 
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For  the  perfect  magnetogas  we  have  described  above, 
these  quantities  are  governed  by  the  Lundqulst  equations 
[8],  which  are  as  follows: 

|H-  +  V.(pu)  =:  0  (1-1-1) 

ll  +  (u.V)u  +  i  Vp  +  i  (B  X  J)  =  0         (1-1.2) 

ll  +  V  K  E  -  0  (1-1-3) 

V-B  =  0  (1-1.^) 

V  X  B  -[XqJ  =  0  (1-1-5) 

E  +  u  X  B  =  0  (1-1.6) 

p  =  p(p)  (1-1-7) 

Equation  (l-l.l),  the  continuity  equation  of  ordinary 
gasdynamics,  is  a  statement  of  the  principle  of  mass 
conservation.   The  next  equation,  a  statement  of  Newton's 
second  law,  is  a  vector  equation  known  as  the  ponderomotive 
equation,  the  last  term  being  the  Lorentz  force.   Equation 
(1-1.3)  is  Faraday's  law  of  induction.   Equation  (l-l.'^) 
which  expresses  the  absence  of  magnetic  charges,  is  not 
independent  of  equation  (1.1-3)-   Since  the  divergence  of  a 
curl  is  identically  zero,  the  divergence  of  equation  (1-1.3) 
is 

^  (V-B)  =  0  .  ,  (1-1.8) 

Thus,  we  can  regard  equation  (l-l.'^-)  as  an  initial  condition 
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rather  than  an  independent  equation  of  motion.   Ampere's 
law,  neglecting  displacement  current,  is  expressed  by 
equation  (I-I.5).   Equation  (1-1.6)  is  Ohm's  law  for  a 
perfect  conductor.   Finally,  equation  (1-1. 7)  is  the  adiabatic 
equation  of  state  for  the  fluid.   Thus,  the  Lundqulst  equa- 
tions Involve  the  coupling  of  Maxwell's  equations  (without 
displacement  current)  with  the  classical  equations  of  gas- 
dynamics.   Although  the  Maxwell  equations  with  displacement 
current  are  Lorentz  invariant,  the  Maxwell  equations  without 
displacement  current  are  Galilean  invariant  [3]-   Since  the 
gasdynamic  equations  are  likewise  Galilean  invariant  (and  we 
are  interested  in  nonrelativistic  motions)  it  is  correct  to 
neglect  the  displacement  current  in  Ampere's  law. 

We  shall  now  eliminate  E  and  J  from  these  equations. 
First,  we  solve  equation  (I-I.5)  for  J  obtaining 

J  =  -f-  (V  X  B)  (1-1-9) 

•^o 

and  then  we  substitute  this  value  into  equation  (1-1.2). 
Similarly,  from  equation  (1-1.6)  we  have 

E  =  B  >:  u  (1-1.10) 

and  we  can  use  this  value  in  equation  (1-1.3)-   These  two 
substitutions  reduce  the  system  of  equations  (l-l.l)  through 
(1-1. 7)  to  the  following  system: 

|^+  V-(pu)  =  0  (1-1-11) 
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|H  +  (u-V)u+i  Vp+i  (B  x[A  (V  X  B)])  =  0      (1-1.12) 

||+Vx(BKu)=0  (1-1.13) 

V-B  =  0  (1-1-1^) 

p  =  p(p)  (1-1.15) 

This  system  of  equations  is  nonlinear.   In  order  to 
linearize  the  system,  we  shall  have  to  restrict  our  discussion 
to  small  disturbances  away  from  the  rest  state,  or,  equiva- 
lently,  to  small  expansion  velocities  of  the  piston.   With 
this  restriction  in  mind,  we  can  linearize  the  equations 
(1-1.11)  through  (1-1.15)  by  setting 

P  =  Po(l  +  P')  (1-1.16) 

u  =  u'  (1-1.17) 

B  =  B  (e,  +  B' )  (1-1.18) 

where  p   and  B  =  Be.,  are  the  constant  undisturbed  values 
of  the  fluid  density  and  the  magnetic  field.   (e'.,  is  a  unit 
vector  which  is  perpendicular  to  the  axis  of  the  piston.)   We 
now  insert  (l-l.lo),  (I-I.17),  and  (I-I.I8)  into  the  nonlinear 
equations,  and  neglect  products  of  primed  quantities.   If  we 
then  drop  the  primes  for  notational  convenience,  the  resulting 
linearized  Lundquist  equations  are: 

|^+  V-u  ^  0  (1-1.19) 


1^  +  a^Vp  +  A^^^  X  (V  X  B)  =  0  (1-1.20) 

ll  +  V  X  (e^  X  u)  =3  0  (1-1-21) 

plus  the  Initial  condition 

V-B  =  0  .  (1-1.22) 

In  these  linearized  equations   a   and  A  are  the  undisturbed 
values  of  the  sound  speed  and  the  Alfven  speeds  respectively. 
They  are  defined  to  be 

(1-1.23) 
(1-1.24) 

If  we  now  introduce  Cartesian  coordinates,  with  the 

z-axis  coincident  with  the  axis  of  the  piston,  none  of  the 

variables  will  exhibit  a  z-dependence .   Thus,  we  can  consider 

the  linearized  Lundquist  equations  in  two  spatial  variables, 

namely  x  and  y.   To  be  definite  we  shall  take  the  x-direction 

to  be  the  direction  determined  by  the  unit  vector  e-,  .   We 

shall  write  u  =  (u,v,w)  and  B=  (B  ,B  ,B  ).   The  equations  for 

X   y   z 

_i     _\ 

the  z -components  of  u  and  B  are 

|^-a2|^=o  (1-1.25) 

s/-i=°  (1-1-26) 


2 
a 

_   dp 
-  dp 

A^ 
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These  two  equations  govern  the  propagation  of  the  Alfven 
wave.   It  Is  easy  to  see  that  equations  (I-I.25)  and  (1-1.26) 
reduce  to  two  wave  equations,  producing  waves  travelling  In 
the  x-dlrectlon  with  speed  A.   We  shall  not  discuss  the 
Alfven  wave  further. 

The  remaining  scalar  equations  of  the  system  (I-I.I9) 
through  (1-1.21)  depend  only  on  the  variables  x,  y,  and  t. 
and  are  thus  the  desired  two-dimensional  formulation.   These 
five  scalar  equations  are  as  follows: 

3u  ^  ^2  a 


St-"  ^'^  (1-1-28) 

|^.a2|£+A2(3-i.^).0  (1-1.29) 


3t^+S7=°     '  (l-l-'°) 

S^-||^°  (1-1-51) 

The  Initial  condition  (1-1.22)  becomes 
Sb    Sb 

It  Is  this  system  of  linear  equations,  depending  on  two 
spatial  variables,  with  which  we  shall  concern  ourselves. 
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1-2.   Plane  Waves 

By  a  plane  wave  we  mean  a  sol-utlon  of  the  equations 
(1-1.27)  through  (1-1.32)  in  which  the  variables  depend 
upon  X,  y,  and  t  only  in  the  combination 

K=Vt-kx-kk  (1-2.1) 

X     y 

In  (1-2.1)  V  is  the  propagation  speed  of  the  wave,  and 

k  =  (k  ,k  )  is  the  unit  vector  in  the  direction  of  propagation 
X  y 

With  the  x-direction  chosen  to  coincide  with  the  direction  of 

B  =  B  t^  ,    and  the  z-direction  chosen  to  coincide  with  the 
o    o  1 

axis  of  the  piston,  the  y-direction  is  determined  by  the 
requirement  that  (x,y,z)  form  an  orthogonal  triad.   We  choose 
1?  to  lie  in  the  xy-plane.   If  we  now  denote  by  9   the  angle 
between  the  two  vectors  k  and  e,  we  have 

k  =  cos  9  (1-2.2) 

X 

k  -  sin  9  (1-2.3) 

y 

To  study  the  plane  wave  solutions  of  the  system  of 

equations  (1-1.27)  through  (1-1.32)  we  make  the  substitutions 

p(x,y,t)  =  p(k),  u(x,y,t)  =  u(k),  v(x,y,t)  =  v(k), 

B  (x,y,t)  =  B^(k),  and  B^^(x,y,t)  =  B  (k),  where  k  is  given 
X  X  N(  y 

by  (1-2.1).   With  these  substitutions  we  obtain  the  following 
system  of  equations: 
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Vp'(k)  -  k  u'(k)  -  k  v'(k)  =  0  (1-2.^) 

Vu'(k)  -  a^k  d'(k)  =  0  (1-2.5) 

Vv'(k)  -  a^kyP'(K)  -  A^(kyB^  -  k^Bp  =   0   (1-2.6) 

VB'(k)  -  k  v'(k)  =  0  (1-2.7) 

X       y 

VB'(k)  +  k  v'(k)  =  0  (1-2.8) 

y     ■'^ 

where  the  primes  denote  differentiation  with  respect  to  the 
variable  k.   This  is  a  homogeneous  system  of  linear  algebraic 
equations  in  the  derivatives  p'(k),  ■u'(k),  v'(k),  B^(k)'  ^-^^ 
B'Tk).   In  order  that  this  system  possess  non-trivial  solu- 

y 

tions  the  determinant  of  the  coefficients  must  vanish.   Using 
equations  (1-2.2)  and  (1-2.3),  this  requirement  gives  the 
dependence  of  the  propagation  speed  V  on  the  angle  9,    namely 

V^  -  (A^  +  a^)V^  +  A^a^cos^e  =  0  (1-2-9) 

Equation  (1-2.9)  has  two  positive  roots,  V^  and  V  , 

I        s 

2 
with  V^  >  V  ,  equality  holding  only  when  A  =  a  and  cos  0=1. 

The  corresponding  waves  are  called  fast  and  slow  waves 

respectively. 

The  dependence  of  V  on  the  angle  0  can  be  illustrated 

graphically  [^].   The  resulting  figure  is  known  as  the 

normal  speeds  locus.   Figure  1  is  the  normal  speeds  locus 

for  the  case  A  >  a.   The  normal  speeds  locus  can  easily  be 

constructed  for  other  cases  [4]. 
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We  shall  now  consider  the  totality  of  plane  waves  which 
can  propagate  in  all  directions  in  the  xy -plane,  and  shall 
suppose  that  the  front  of  each  wave  passes  through  the  origin 
at  time  t  =  0.   Then,  at  time  t  =  t  ,  the  various  wave  fronts 
form  two  families  of  straight  lines.   The  envelope  of  these 
straight  lines  (characteristics)  is  the  t  =  t   cross-section 
of  the  characteristic  conoid.   The  characteristic  conoid 
gives  the  range  of  influence  of  a  disturbance  at  the  origin 
at  time  t  =  0.   We  can  derive  the  parametric  equations  for 
the  curves  which  make  up  the  characteristic  conoid  in  a 
straightforward  manner  [4]. 

Suppose  that  an  initial  disturbance  is  on  a  plane  X    (6) 
and  proceeds  in  the  normal  direction  given  by  the  angle  9 
(see  figure  2).   After  a  time  t  the  disturbance  has  moved 
to  a  plane  X, (0)  having  traversed  a  distance  V^(9)t  where 
n  -   s  or  f,  depending  on  whether  the  disturbance  is  a  slow 
or  fast  wave.   We  have  assumed  that  the  plane  X  (0)  passes 
through  the  origin,  and  thus  can  write  the  equation  of  this 
plane  as 

X  cos  0  +  y  sin  0=0.  (1-2.10) 

The  equation  of  the  plane  X,  (0)  is  then 

X  cos  0  +  y  sin  0  =  ±  V  (0)t  (1-2. 11) 

where  the  plus  sign  corresponds  to  the  direction  of  propagation 
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shown  In  figure  2,  while  the  mln-us  sign  corresponds  to  the 
opposite  direction  of  propagation. 

The  envelope  of  the  planes  X,  (9)  Is  given  by  the 
values  of  (x/t)  and  (y/t)  which  satisfy  the  equation 
(1-2.11)  and  Its  derivative  with  respect  to  Q,    namely, 

dV  (0) 
-  X  sin  0  +  y  cos  e  =  t  ( — ^ — )t  (1-2.12) 


These  values  are 


,dv^(e) 


r  =  t   [V^(0)cos  0  -  (  "^  ')  sin  0]       (1-2.13) 


dV  (0) 

I  =  t    [V^(0)sln  0  +  (   gg   )  cos  0]  (1-2.14) 

We  can  determine  explicit  expressions  for  V„(0)  and  V  (0) 

I  s 

by  solving  equation  (1-2.9).   These  expressions  are 

V^(0)  =  A[|  (1+s)  +  y(0)]i  (1-2.15) 

V3(0)  =  A[|  (1+s)  -  y(0)]^  (1-2.16) 


where 


and 


,A  2 


s  =  (f)  (1-2.17) 


Y(0)  =  [(l+s)2  -  4s  cos^0]2  (1-2.18) 

Using  equations  (1-2.15)  through  (1-2.18)  we  find  that 
the  fast  wave  envelope  is  given  by  the  parametric  equations 
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f  =  t   Acos    e[a{9)    -^jl^]  (1-2.19) 

l^t  A  sm  G[a{e)  +  ICoMo)  ^  ('-"-'^^ 


where 


a{e)    =  -^-^—  (1-2.21) 

The  slow  wave  envelope  is  given  by  the  parametric  equations 

|  =  ±  ACOS  8[p(e)  +f^|^]  (1-2.22) 

2=t  Aslne[p(9)  -f^fl^]  (1-2.25) 

where 

V  (0) 
P(0)  =  -^  (1-2.24) 

In  figure  3  we  have  sketched  the  characteristic  conoid, 
and  in  figure  4  we  have  sketched  a  t  =  constant  cross-section 
of  the  conoid  [4].   The  outer  branch,  ABCD,  is  the  fast  wave 
envelope.   The  slow  wave  envelope  is  composed  of  the  two 
cusped  curves,  cFaE  and  GIKH,  and  the  straight  line  segment 
Gc.   The  distance  Oc  is  given  by 


^^  (1-2.25) 


/  A  +  a 
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1-5.   Representation  of  a  Conical  Field  by  Plane  Waves 

We  are  going  to  look  for  a  solution  of  equations 
(1-1.27)  through  (1-1.32)  which  is  self-similar,  or  conical. 
By  this  we  mean  that  the  variables  will  depend  only  on  x/t 
and  y/t,  or,  put  another  way,  the  variables  will  be  homo- 
geneous functions  of  x,  y,  and  t  of  degree  zero.   This  means 
that  snapshots  of  the  flow  at  successive  instants  would  be 
one  picture,  being  enlarged  uniformly  at  a  uniform  rate.   We 
shall  construct  such  a  solution  by  superposing  plane  waves 
that  are  conical. 

The  solution  to  equations  (I-I.27)  through  (1-1.32)  can 
be  written  in  terms  of  a  single  function  if {x,y,t),    [5]-   In 
particular,  equation  (1-1.32)  implies  the  existence  of  a 
function  ^(x,y,t)  with 

B      =   +    I  (1-3.2) 

y       ^X 

where  the  subscripts  appended  to  the  function  ^  denote 
partial  differentiation.   Gardner  [5]  realized  that  the  choice 

J(x,y,t)  =  ^y^(x,y,t)  (1-3.3) 

enables  all  the  variables  to  be  expressed  in  terms  of  simple 
combinations  of  third  derivatives  of  the  function  lA.  Refer- 
ring back  to  equations  (1-1. 27)  through  (I-I.32)  we  find 
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that  (1-3.1)  through  (1-3.3)  Imply  that 


Si 


V  .  -  ^y,,  (1-3.6) 


B  =  -  -^^  ,  (1-3.7) 

X      yyt 


B  =  +  ^  ,  (1-3.8) 

y      xyt 


Equation  (1-1.27)  then  yields  the  following  equation  for  the 
function  i'{x,Y ,t)  : 

St  St   Sx    hy 

^   ^2,2  S^   /S^V'  ^  S^^N    ^  fi    -^  o\ 

+  a  A  — ^  ( — 2  +  — 2^  "^  (1-3.9) 

Sx    Sx    Sy 


Since  we  are  trying  to  construct  a  conical  flow,  we  want 
the  third  derivatives  of  i/   to  be  homogeneous  functions  of 
X,  y,  and  t  of  degree  zero.   Thus  if   itself  must  be  a  homo- 
geneous function  of  x,  y,  and  t  of  degree  three.   Referring 
to  formula  (1-2.1)  for  the  form  of  a  plane,  we  write 

z 
i'{x,y,t)    =  j  (t  -  A(s)x  -  [x(s)y)\(s)ds    (1-3-10) 

z 

o 
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We  write  1^   in  the  form  of  an  integral  so  that  differentia- 
tion will  be  particularly  simple.   The  point  z^  is  fixed. 
The  point  z  is  a  function  of  x/t   and  y/t  such  that  the 
integrand  vanishes  at  the  point  z,    i.e., 

t  -  A(z)x  -  n(z)y  =  0  (1-3.11) 

Differentiating  this  relation  gives 

^:^:^=  1  :  -A  :  -n  (1-3-12) 

dt   ox   oy 

With  the  aid  of  these  results  we  can  compute  the  various 
derivatives  of  the  integral  representing  the  function  i/(x,y,t) 
The  expression  (1-5.10)  will  satisfy  the  equation  (1-3-9) 
provided  that 

1  -  (a^  + A^)(a2  +  h2)  +a^A^A2(A2+n2)  =  o    (1-3-13) 

The  significance  of  A  and  M-  can  be  seen  by  examining 
(1-3-11)  and  (1-5-13)-   First,  we  see  that  (1-3-11)  is  the 
equation  of  a  straight  line  in  (tt?  p) -space,  and  that  l/A 
and  I/m-  are  the  Intercepts  of  this  line  with  the  x/t  and 
y/t  axes  respectively.   Equation  (1-3.15)  Implies  that  this 
straight  line  is  a  tangent  through  the  point  P  =  (y,    y)    to 
the  characteristic  conoid  of  figure  5-   In  general,  through 
a  given  point  P,  there  are  four  such  tangents.   (Upon 
eliminating  i-l  from  (1-5. 15)  with  the  aid  of  (1-5. 11),  we 
obtain  a  quartic  equation  for  A.)   The  nature  of  the  four 
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tangents  depends  upon  the  location  of  the  point  P,  [9,10]. 

For  a  point  P-,  outside  the  outer  branch  of  the  char- 
acteristic conoid  the  Intercepts  of  these  four  tangents  are 
real,  and  the  corresponding  waves  are  two  fast  waves  and 
two  slow  waves.   For  a  point  Pp  between  the  Inner  and  outer 
branches  of  the  characteristic  conoid,  two  of  the  tangents 
have  real  Intercepts,  and  the  corresponding  waves  are  slow 
waves.   The  other  two  tangents  from  such  a  point  have  complex 
Intercepts.   These  two  cases  are  Illustrated  In  figure  5- 
Finally,  for  a  point  P^  Inside  the  Inner  branch  of  the  char- 
acteristic conoid,  the  Intercepts  of  all  four  tangents  are 
real,  and  the  resulting  waves  are  four  slow  waves.   This 
situation  Is  shown  In  figure  6.   We  shall  be  concerned  with 
the  region  between  the  Inner  and  outer  branches  of  the 
characteristic  conoid.   We  shall  refer  to  the  solution  corres- 
ponding to  the  two  tangents  with  real  Intercepts  as  the 
hyperbolic  wave,  and  the  solution  corresponding  to  the  two 
tangents  with  complex  Intercepts  will  be  called  the  elliptic 
solution. 

In  accordance  with  (1-3. 10)  we  shall  write  our  solution 

(formally)  In  the  form 

z 
p  =  j   P(s)ds  (1-3-1^) 

^o 
z 

u  =  I   U(s)ds  (1-3.15) 

^o 
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z 

V 

z 


z 

=  [  V(s)ds  (1-3.16) 


o 


\ 


Zi 

=  f  B3_(s)ds  (1-3.17) 


z 

o 


z 

By  -  [   B2(s)ds  (1-3.18) 

z 

o 

In  the  hyperbolic  regime, z  will  be  real,  and  these  Integrals 
will  be  real.   In  the  elliptic  solution,  z  will  be  complex 
and  these  integrals  will  be  complex  contour  integrals  of 
analytic  functions.   We  shall  impose  further  conditions  on 
the  integrands,  the  paths  of  integration,  and  choose  the 
point  z  .   We  can  differentiate  these  formal  expressions, 
and  use  equations  (I-I.27)  through  (I-I.32)  to  obtain  rela- 
tions between  P,  U,  B-,  ,  Bp ,  and  V   (V  is  singled  out  for 
convenience  only.)   If  we  use  the  relations  contained  in 
(1-3.12)  we  find  that 


P(s)  =   g  ^    [1  -  A^(a2(s)  +  ^l2(s))]V(s)  (1-3.19) 
a  \i{s) 

U(s)  =   ^  [1  -  a2(a2(s)  +  li^{s))]Y{s)         (1-3.20) 


B^(s)  =  n(s)V(s)  (1-3.21) 

B2(s)  =  -  A(s)V(s)  (1-3.22) 
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1-4.   Jump  Relations 

For  future  reference  we  shall  discuss  in  this  sub- 
section the  propagation  of  discontinuities  in  a  solution 
of  the  equations  (I-I.27)  through  (1-1-32).   Let  us  suppose 
that  we  have  a  solution  of  these  equations  which  tends  to 
different  limits  as  the  point  (x,y,t)  tends  to  the  two 
different  sides  of  a  surface  S  in  (x,y,t)  space,  so  that 

the  quantities  p,  u,  v,  B  ,  and  B  have  Jump  discontinuities 

^      y 

In  such  circumstances  the  differential  equations  (I-I.27) 
through  (1-1. 32)  break  down;  in  their  place  we  have  the 
jump  relations  given  below.   These  jump  relations  can  be 
obtained  in  a  variety  of  ways,  for  instance,  by  considering 
weak  solutions  of  the  equations  under  discussion.   We  denote 
the  components  of  the  normal  vector  to  S  by  n  ,  n  ,  and  n,  . 
The  desired  jump  relations  are  as  follows  [l]: 

n^[p]  +  n^[u]  +  ny[v]  =  0  (l-^-l) 

n^[u]  +  a^n^[p]  =  0  (1-4-2) 

n^[^]    +   a^ny[p]  +  A^(ny[B^]  -  n^[By])  =  0  (1-4.5) 

n^[B^]  +  ny[v]  =  0  (1-4.4) 

n^[By]  -  n^[v]  -  0  (1-4-5) 

n[B]+n[B]=0  (1-4.6) 

X     X-'  y     y 

where  [Q]  denotes  the  jump  in  the  quantity  Q  across  the 
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surface  S.   We  can  eliminate  the  jumps  In  the  five  variables 

from  these  six  equations,  obtaining  an  equation  relating 

the  components  of  the  normal  vector  to  the  surface  S.   This 
equation  Is 

nj-  (A^  +  a2)n2(n2  4-n2)  +  AVn^(n^  +  n2)  =  o     (1-4.7) 

Equation  (l-'^^.T)  Implies  that  the  surface  S  Is  a  character- 
istic surface  of  the  system  of  equations  (1-1.27)  through 
(1-102). 


1-5-   Boundary  Conditions  on  the  Piston 

In  this  subsection  we  shall  discuss  the  form  of  the 
linearized  boundary  conditions  on  the  piston.   We  shall  use 
a  slightly  different  form  of  the  boundary  conditions  than 
those  obtained  here  when  we  discuss  in  the  last  section  pistons 
which  are  nearly  parallel  to  the  undisturbed  magnetic  field 
Bq.   All  vectors  in  this  subsection  are  two-dimensional 


vectors  in  {tt,   ^) -space 


We  shall  impose  only  one  restriction  on  the  shape  of  the 
piston,  namely  that  (n-e,)  =  0  at  only  a  finite  number  of 
Isolated  points.   Here  n  is  the  unit  normal  to  the  surface 
of  the  piston,  and  e,  is  a  unit  vector  in  the  ^  direction. 
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This  assumption  permits  us  to  cast  the  boundary  conditions 
into  a  useful  form.   We  shall  denote  by  U   the  velocity 
with  which  the  piston  is  expanding. 

The  first  boundary  condition  to  be  applied  at  the  sur- 
face of  the  piston  is  that  no  fluid  flows  through  that 
surface ,  i.e., 

u-n  =  U  -n  (1-5.1) 

We  shall  denote  the  magnetic  field  on  the  surface  of  the 
piston  by  B  ,  and  tlie  corresponding  electric  field  by  E  . 
We  assume  that  on  the  surface  of  the  piston 

E  +  U  >^  B  =0  (1-5-2) 

p    p    p  \      ^      I 

This  assumption  is  equivalent  to  the  requirement  that  the 
surface  of  the  piston  is  a  perfect  conductor.   Using  Ohm's 
law  for  the  perfectly  conducting  fluid,  (1-1.6),  we  can 
express  the  condition  that  the  tangential  component  of  the 
electric  field  is  continuous  across  the  surface  of  the 
piston  by  writing 

n  >^  (B  X  u)  =  n  ^  (B"  ^  UL )  (1-5-3) 

We  shall  now  linearize  equation  (1-5.3)-   In  analogy 
with  (1-1. 18)  we  write 

B   =  B  (e.  +  b' )  (1-5-^) 

p     o''  1    p^  y      -^      I 
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If  we  denote  by  dS  the  element  of  surface  area  on  the 
piston,  we  know  that 


Dt 


/ 


/• 


B-dS 


■SB 


(^+  V  X  (B  Xu)).dS        (1-5.5) 

s(t) 

The  nonlinear  equation  (1-1.13)  Implies  that  the  Integrand 
In  the  right  hand  Integral  Is  zero,  so  that  on  any  piece  of 
the  piston,  S(t),  we  must  have 

^•dS'  =  constant  (I-5.6) 

s'(t) 

For  a  conical  flow,  B  remains  constant  while  S(t)  changes 
with  time.   Thus,  the  constant  In  (I-5.6)  must  be  zero. 
This  can  be  achieved  by  taking  B  to  be  tangential  to  the 
surface  of  the  piston.   Thus,  one  might  suspect  that  the 
linearization  (1-5.4)  would  be  justified  only  for  a  piston 
whose  surface  Is  nearly  parallel  to  the  undisturbed  magnetic 
field.   This  conclusion  Is  borne  out  by  our  subsequent 
analysis.   It  should  be  pointed  out,  however,  that  the 
a  priori  conclusion  Is  based  upon  the  nonlinear  equations; 
this  conclusion  cannot  be  drawn  from  the  linearized  system. 
A  priori,  the  linearized  equations  and  boundary  conditions 
are  consistent. 

Using  equations  (I-I.I8)  and  (1-5-4)  in  (I-5.3)  we 
obtain  the  linearized  boundary  condition 
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n  X  (e^  ?^  u)  =  n  X  (e^  X  U  ^ )  (1-5.?) 


Upon  expanding  the  vector  triple  products  we  find  that 

(n.v)e^  -  (n.e^)u  =  (n-Up)e^  "  (n-e^)Up  (1-5-8) 

We  can  rearrange  this  equation  to  obtain 

(n.e^)(u  -  ifp)  =  0  (1-5-9) 

With  our  assumption  on  the  vanishing  of  (n-e. ),  and  the 
desired  continuity  of  u,  we  conclude  that 

u  =:  Up  (1-5-10) 

on  the  surface  of  the  piston.   Thus,  the  requirements  that 
no  fluid  flow  through  the  surface  of  the  piston,  and  that 
the  tangential  component  of  the  electric  field  Is  continuous 
across  that  surface  lead  to  the  two  boundary  conditions 

u  =  Up-e^  =   U3_  (1-5-11) 

V  =  Up-^2  =  ^2  (1-5-12) 

at  the  surface  of  the  piston.   (ep  is  a  unit  vector  in  the 
y-dlrection . ) 

Finally,  we  must  require  that  the  normal  component  of 
the  magnetic  field  is  continuous  across  the  surface  of  the 
piston.   In  terms  of  the  perturbations  of  the  magnetic  field 
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on  the  piston  and  In  the  fluid  we  require  that 

B-n  =  B  -rT  (1-5.13) 

This  relation  provides  the  definition  of  B-n  on  the  surface 
of  the  piston.   We  can  then  solve  an  interior  problem  to 
obtain  the  magnetic  field  throughout  the  piston. 

In  the  last  section,  when  we  consider  pistons  which  are 
everywhere  nearly  parallel  to  the  undisturbed  magnetic  field, 
we  shall  consider  a  slightly  different  form  of  these 
boundary  conditions. 


1-6.   Uniqueness  of  the  Solution 

We  shall  model  our  uniqueness  theorem  upon  that  given 
by  Gardner  [5].   We  suppose  that  p,  u,  v,  B  ,  and  B  are 

^      y 

differentiable  in  the  open  region  R  of  (x,y ,t ) -space  which 
is  outside  the  piston  and  does  not  include  certain  surfaces, 
and  that  these  quantities  satisfy  the  equations  (I-I.27) 
through  (1-1.32)  in  R.   The  surfaces  include  the  surface  of 
the  piston,  the  characteristic  surfaces  ABCD,  cFaE,  GIKH, 
and  that  portion  of  the  line  segment  Gc  which  lies  outside 
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of  the  piston,  and  certain  line  segments  across  which 
the  variables  are  continuous,  but  their  various  first 
derivatives  are  not  (see  figure  7)-   We  assume  that  at  a 
characteristic  surface  the  jump  relations  (l-^.l)  through 
(1-4.6)  hold.   Finally,  we  assume  that 


p=u=v=B^=By=0  (1-6.1) 


for  t  >  t   outside  some  bounded  region  of  the  xy-plane. 

This  is  a  statement  of  the  fact  that  magnet ohydrodynamic 

waves  propagate  with  finite  velocity. 

To  show  that  such  a  solution  is  unique,  we  suppose  that 

there  could  be  two  solutions,  p,  ,  u-,  ,  v-,  ,  B  -.  ,  B  -,  and 

p„ ,  u„ ,  v^ ,  B  o J  B  o •   We  shall  consider  the  difference  of 
^2         2 '   2    x2 '   y2 

these  two  solutions:   p  =  p -.  -  Pp ,  u  =  u-,  -  Up ,  v  =  v-,  -  Vp , 

B  =  B  -,  -  B  oj  and  B  =  B  -,  -  B  ^  •   We  shall  now  show 
X     xl     x2'       y      yl     y2 

that  this  solution  is  identically  zero. 

The  proof  employs  the  abc  or  energy  method.   For  this 

method  to  succeed,  the  various  variables  must  be  square 

integrable.   We  multiply  (1-1.27)  by  a  p,  (1-1.28)  by  u, 

(1-1.29)  by  V,  (1-1.30)  by  A^B   and  (l-l.^l)  by  A^B  .   We 

X  y 

then  add  the  six  resulting  equations  to  obtain  a  relation 
which  may  be  written  in  the  form 

1  ^  (a2p2  +^;2  ^^2^  ^2^2  ^  ^2^2^  ^  J_  (a2p^  _  ^^vB^) 

+  -I-  (a^pv  +  A^vB^)  =  0  (1-6.2) 


-^   ^ct  pv   Ih  H   VID^, 
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We  now  Integrate  (1-6.2)  over  all  (x,y,t)  space,  exclusive 
of  the  piston,  with  t  <  t^.  On  the  surface  of  the  piston, 
the  solution  will  satisfy  the  homogeneous  boundary  conditions 

u  =  V  =  B-n  =  0  (1-6.3) 

This,  combined  with  the  square  integrability  of  the  solution, 

insures  that  the  boundary  terms  arising  from  the  integration 

at  the  surface  of  the  piston  will  vanish. 

Moreover,  the  boundary  terms  arising  at  a  characteristic 

surface  S  likewise  contribute  nothing.   If  we  denote  the 

values  of  the  variables  on  one  side  of  S  by  p ,  u ,  v ,  B^ ,  and 

B  ,  and  the  values  on  the  other  side  of  S  by  p ,  u,  v ,  B^, 
y  ^ 

and  B  ,  then  the  integrand  of  the  surface  integral  arising 
at  S  is 

1      'p    '^       _P    P    P— P    P— ?         P  —    2 

i  n,  (a  p  +  u-  +  V  +  A  B^  +  A  B  )  +  n^(a''pu  -  A  vB^  ) 
d      z  X     y     X  y 

+  n  (a  pv  +  A  vB  )  -  ^  n,  (a^p^  +  u  +  v  +  A  B^  +  A^'b^J 
y  X    £1   13  X      y 

-  n^(a^pu  -A^vBy)  -  ny(a^pv  +  A^vB^)  (1-6.4) 

We  can  regroup  equation  (1-6.4)  to  obtain  the  expression 
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|(p  +  p)  f  n^[p]  +n^[u]  +ny[v]  j  +  |(u  +u)|n^[u]  +a^n^[p]| 
+  |(v+v)|n^[v]  +  a^ny[p]+A^ny[B^]  -  A^n^[By]  j    (1-6. 5) 

Referring  to  the  jump  relations  (1-4.1)  through  (1-4-.6)  we 
see  that  (1-6. 5)  is  identically  zero. 

Finally,  the  integrated  terms  arising  at  a  line  where 
the  solution  is  continuous  and  its  derivatives  are  not  also 
vanish.   (Such  a  line  is  also  a  characteristic,  but  not 
part  of  the  characteristic  conoid.)   Therefore,  upon  inte- 
grating (1-6.2)  we  are  left  with  the  following  integral: 


1_ 
2 


(a^p  +  u^  +  v^  +  A^B^  +  A^B^)dxdy  =0   (1-6.6) 

X     y 


t=t^ 

Equation  (1-6.6)  implies  that 

p=u=v=B^-By=0  (1-6.7) 

which  is  what  we  wanted  to  show. 
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2.   The  Elliptic  Solution 

In  the  first  section  we  discussed  the  nature  of  the 
tangents  to  a  t  =  constant  section  of  the  characteristic 
conoid.   In  particular,  for  a  point  P  between  the  slow  and 
fast  wave  envelopes,  we  saw  that  two  of  the  four  tangents 
to  the  characteristic  conoid  had  complex  -^   and  ^-intercepts. 
It  is  these  characteristics  that  are  employed  in  the  con- 
struction of  the  elliptic  solution. 

In  the  first  subsection  we  obtain  an  explicit  repre- 
sentation of  the  intercepts  of  these  tangents,  j-   and  — ,  in 
terms  of  a  complex  parameter  z.   The  relation  between  z  and 
Y   and  y   provides  a  non-conf  ormal  mapping  of  the  (^,  |^) -plane 
into  the  complex  z-plane.   Following  Gardner  [5]^  we  devote 
the  next  subsection  to  a  detailed  examination  of  this 
mapping.   By  employing  this  mapping,  we  formulate  our  problem 
as  a  boundary  value  problem  for  an  analytic  function  in  the 
z-plane . 

In  the  final  subsection  of  this  chapter  we  discuss 
this  formulation  of  the  problem,  and  its  solution.   The 
solution  is  completely  determined  up  to  one  analytic  function. 
Certain  general  restrictions  are  placed  on  this  function,  but 
its  complete  determination  depends  upon  the  boundary  condi- 
tions imposed  on  the  surface  of  the  piston,  and  upon  the 
shape  of  the  piston  cross-section. 
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2-1.   The  Form  of  A  and  U- 

In  the  first  section  we  wrote  our  solution  In  terms 
of  various  third  derivatives  of  a  homogeneous  function  of 
X,  y,  and  t  of  degree  three,  namely 


■^(x 


z 

,y,t)  =  [  (t  -A(s)x  -n(s)y)\(s)ds     (2-1.1) 


z 

o 


The  formal  solution  was  then  given  by  equations  (1-5-'^) 
through  (1-3.8).   We  concluded  at  that  point  that  A  and  \i 
must  also  satisfy  the  two  relations 

t  -  A(z)x  -  n(z)y  -  0  (2-1.2) 


1  -  (a^  +A^)(A^(z)  +n^(z))  +  a^A^A^(z)(A^(z)     (2-1.3) 

+  M-^(z))  =  0 


If  we  Insert  a  specific  form  for  i^  (z)  Into  (2-1.3) 

2 
we  obtain  a  quadratic  equation  for  A  (z).   If,  on  the  other 

2 
hand,  we  substitute  a  specific  form  for  A  (z)  Into  (2-1.3) 

2 
we  obtain  a  linear  equation  for  li-  (z).   We  follow  this 

latter  course  of  action  and  write  [5] 

A(z)  =  1  (2-1.4) 

Geometrically,  z  Is  the  complex  (ir) -Intercept  of  a  tangent 
to  the  characteristic  conoid.   Using  (2-1. •4-)  In  (2-1.3)  we 
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obtain 


^(z)    =    (^jifi^lApiTIZi  (2_i.5) 

In  order  that  H.(z)  be  a  single -valued  function  of  the 
complex  variable  z,  we  must  make  certain  cuts  In  the 
z-plane.   We  do  this  by  making  four  cuts  on  the  real  axis, 
from  -A  to  -oo,  from  -c  to  -a,    from  c  to  a,  and  from  A  to  oo 
(see  figure  8  ).   To  choose  a  definite  branch  of  the  square 
root  In  (2-1.5)  we  write 


,(,)  .iVd-^Vm-Z/a")  (2-1.6) 

^  '       1  -tT/q. 


and  choose  this  square  root  to  be  +1  for  z  =  0. 

The  solutions  which  result  from  consideration  of  the 
two  complex  characteristics  are  complex  themselves.   We  want 
to  form  a  combination  of  these  two  solutions  which  is  real. 
Since  the  two  complex  solutions  are  complex  conjugates  of 
one  another,  we  shall  consider  half  their  difference,  or 
equlvalently,  the  Imaginary  part  of  either  solution.   We 
can  accomplish  this  by  modifying  formulae  (2-1.1).   We  take 
as  our  definition  of  TA(x,y,t)  for  the  elliptic  solution 
the  formula 

z 

•^(x,y,t)  =  Im  J  (t  -  A(s)x  -  n(s)y)\(s)ds   (2-1.7) 

z 
o 
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In  this  formula  we  shall  take  h(s)  to  be  an  analytic 
function.   The  Integral  will  be  Interpreted  as  a  contour 
integral  along  an  appropriately  specified  contour. 

Equation  (2-1.2)  gives  the  relation  between  z  and 
jr   and  -p,  namely 


(f )  ^  (^) 


K 


2    2       ^t' 
z     -  c 


(2-1.8) 


By  considering  the  real  and  imaginary  parts  of  this  equation 
we  find  that 


i-m 


> 


Im(z) 


(z  -  A  )  (z  -  a  ) 

2    2 
z  -  c 


Im 


> 


I  =  «-(-)  -  & 


■Rey(-'-^^)(V^S](Z) 


z  -  c 


(2-1.9) 


(2-1.10) 


Equation  (2-1.8)  provides  a  non-conf ormal  mapping  of  the 
[^,   1^) -plane  into  the  complex  z-plane.   We  shall  consider 
this  mapping  in  depth  In  the  next  subsection. 
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'      z  -  c 


2-2.   The  Mapping  (f)  +  (C^y^z  -A  Mz  -a  ;  ^Z) 


As  we  pointed  out  In  the  last  subsection,  this  Is  a 
non-conf  ormal  mapping  of  the  (r-,  #-) -plane  Into  the  z -plane. 
This  mapping  possesses  some  very  special  properties  which 
will  aid  greatly  In  the  construction  of  solutions  to  our 
problem.   In  this  subsection  we  shall  prove  that  this 
mapping  possesses  the  following  five  properties: 

(l)  The  mapping  (2-1.8)  of  the  silt  z-plane  of 

figure  8  Is  one-one  and  dlfferentlable  Into  a 
region  R  of  the  (x-,  ^) -plane . 

(11 )  Points  In  R  for  which  x-  =  0  are  Invariant  under 
this  mapping. 

(ill)  Near  the  origin,  the  mapping  Is  approximately  the 
Identity  mapping;  more  precisely 

2  =  (|)  +  l(f)[l  +  0(|z|^]  (2-2.1) 

for  small  \z\ 

(iv)  The  boundary  curves  of  R  determine  characteristic 
surfaces . 

(v)  A  point  on  the  boundary  of  R  Is  mapped  Into  the 
point  P  on  the  ^--axls  where  the  tangent  from 
that  point  on  the  boundary  of  R  Intersects  the 

X 

•ir-axis . 
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Proof  of  (1):   We  take 

z  =  ^  +  ni  (2-2.2) 

and  write 

A(z)  ==  \{i,r])   +    lA.(|,Ti)  (2-2.3) 

H(z)  =  H^(4,ri)  +  l\i^{i,r\)  {2-2A) 

Recalling  the  representations  (2-1.4)  and  (2-1.5)  for  A 
and  p.,  we  can  write  equations  (2-1.9)  and  (2-1.10)  In  the 
form 

M-. 

t  =  A  ^L.  -  \.[i  (2-2.5) 

r  1    1  r 

A. 

f  -  -  A  ^.  -  A.^  (2-2.6) 

r  1    1  r 

We  shall  now  verify  that  the  denominators  In  these  expres- 
sions to  not  vanish  In  the  silt  z-plane  of  figure  8. 
We  define 

D(|,ri)  =  A^(^,T])^x.(|,Ti)  -  A.(^,r|)n^(l,Ti)   (2-2.7) 
which  can  be  written  In  the  more  convenient  form 

D(^.n)  =    -(4   +  ^^l)^"^[^]  (2-2.8) 
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or,   ivith  the   aid   of    (2-1.4)   and    (2-1.5)   we  may  write 


D{i,r]) 


-if)\u.i^) 


Im' 


2        2 
z     -  c 


(z2-a2)(z2_A2) 


(2-2.9) 


Now  |m-(z)|   Is  everywhere  positive  In  the  silt  z-plane.   We 
will  show  that  the  other  factor  In  (2-2.9)  never  vanishes, 
by  applying  the  maximum  principle  for  harmonic  functions. 
We  excise  small  regions  near  z  =  t  a  and  z  =  t  A,  and 
consider  the  harmonic  function  In  the  resulting  domain. 
Near  the  points  z  =  +  a  and  z  =  ±  A  we  have 


Im' 


2 


/^     2n  /  2    7^ 
(z   -  a  )(z   -  A  ) 


<  0 


(2-2.10) 


while  on  the  silts  we  have 


Im 


2     2 
z   -  c 


(z^ 


2^  /  2 
a  )(z 


a') 


=   0 


(2-2. ii; 


Since  an  harmonic  function  can  assume  Its  maximum  only  on 
the  boundary  of  a  domain,  we  conclude  that  (2-2.10)  holds 

throughout  the  silt  z-plane.   Coupled  with  the  fact  that 

2 
|lJ.(z)|   Is  positive,  we  have 


D(^,Ti)  <  0 


(2-2.12) 


In  the  silt  z-plane  of  figure  8.   Thus,  (2-1.8)  defines  a 
dlfferentlable  mapping  of  the  z-plane  Into  the  (^,  |^) -plane. 


-  36  - 


We  will  now  show  that  the  mapping  is  one -one. 

To  show  that  the  mapping  Is  one -one  we  must  show  that 
for  given  ^  and  ~   there  is  not  more  than  one  value  of  z 
such  that  (2-2.5)  and  (2-2.6)  are  satisfied.   The  possible 
values  of  z  are  the  four  roots  of  the  quartic  equation 


[^t(z)(f)]^  =  [1  -  A(z)(|)]2  (2-2.13) 


This  equation  always  has  two  real  roots,  one  positive  and 

2    2    2 
one  negative,  such  that  c  <  z  <  a  .   The  existence  of  the 

negative  root  is  assured  by  the  fact  that  the  family  of 

lines 


(^)  t  (^)V^^'  -^.')(%-^')  (|)  =  z     (2-2.1^^) 
'       z   -  c 

(-a  <  z  <  -c) 

covers  the  whole  (:p,  ^) -plane  (the  j-   Intercept  of  such  a 
line  can  have  all  values  from  -00  to  +00)  .   a  similar  observa- 
tion with  c  <  z  <  a  guarantees  the  existence  of  the  positive 
real  root.   Thus,  (2-2.13)  can  have  only  two  complex  roots 
which  are  complex  conjugates.   By  equations  (2-2,6)  and 
(2-2.12)  we  see  that 

sgn[lm(z)l   =    sgn[|  I  (2-2.15) 

Only   one   of   the    two  complex-conjugate   roots   of    (2-2.13)   will 
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satisfy  equation  (2-2.15).   This  shows  that  the  mapping 
(2-1.18)  Is  one-one,  except  for  the  regions  -A  <  z  <  -a., 
-c  <  z  <  c,  and  a  <  z  <  A.   In  these  cases,  \i{z)    is  pure 
Imaginary  and  A(z)  Is  real.   From  equations  (2-1.9)  and 
(2-1.10)  we  conclude  that 


X 

^  =  t 


(2-2.16) 


which  Is  a  unique  value  of  z 


Proof  of  (11) :   This  result  v/as  proved  above,  but  can 
be  seen  directly  from  equation  (2-1.8). 

Proof  of  (ill) :   A  straightforward  power  series  expan- 
sion yields 


^T?7  4        6 

z    .     z 


'        1  -  z 


/c^  2a^A^   2a-A^c^ 

+  0(|z|^)   (|z|  <  c)       (2-2.17) 

Using  this  result  In  (2-1.6)  then  gives  the  desired  relation, 
namely  (2-2.1). 

Proof  of  (Iv)  :   With  z  =  ^  +  Ir) ,  we  consider  |  <  -A, 
-a<^<  -c,  c<^<a,  orA<^,  and  assume  that  n  Is  small, 
so  that  we  will  be  near  the  boundary  of  the  silt  z-plane. 
With  I  In  these  Intervals  and  t\   =   0 ,    A(z)  and  l^(z)  are  real. 
Expanding  In  powers  of  q  we  find  that  the  real  and  Imaginary 
parts  of  (2-1.8)  are 
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'>(e)(|)  +  V.{1.){1)    =  1  (2-2.18) 

The  sign  of  \.i   in  these  eqiiatlons  will  depend  on  whether 
r\   — >  0  through  positive  or  negative  values.   Equations 
(2-2.18)  and  (2-2.19)  determine  the  image  in  the  (x-,  |^) -plane 
of  the  bounding  slits  in  the  z-plane.   We  will  write  these 
equations  as 

M^)(|)  +  \^{i){l)    =  1  (2-2.20) 

A'(^)(|  +  H'(l)(f)  =  0  (2-2.21) 

where  the  primes  in  (2-2.21)  denote  differentiation  with 
respect  to  ^.  The  determinant  of  these  two  equations  is 
not  zero.  Thus,  the  point  (x,y,t)  is  on  the  envelope  of 
the  characteristic  planes  given  by 

A(^)x  +  ^L(^)y  =  t  (2-2.22) 

To  see  that  the  determinant  of  equations  (2-2. 20-)  and  (2-2.21) 
does  not  vanish,  we  write  the  determinant  in  the  form 


A,.  -A.,=  a2|,(,)  ^A^^f(,)2^      jg_2^2j) 
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and  compute  directly  that  ^  [  (i^/A)  ]  vanishes  only  for 

I  =  0  or  ^  =  t  V^2c",  and  these  points  are  not  on  the  bounding 

silts . 

Proof  of  (v) :   The  line  given  by 

>^(e)(|)  +  ^^(^)(|)  =  1  (2-2.24) 

Is  tangent  to  the  characteristic  conoid  at  the  point 

(|)  =  (f)    and    (|)  =|)  (2-2.25) 

and  the  Intercept  of  this  line  on  the  (:^)-axls  Is  at 

Concisely  stated,  under  the  mapping  (2-1.8)  the  right 
half  of  the  fast  wave  envelope,  BAD  (see  figure  4),  maps 
onto  the  silt  from  A  to  c»  (see  figure  8).   Similarly,  the 
left  half  of  the  fast  wave  envelope,  BCD,  maps  onto  the 
silt  from  -00  to  -A.    The  two  cusped  portions  of  the  slow 
wave  envelope  map  onto  the  two  silts  from  -a  to  -c  and 
c  to  a.   Finally,  distances  along  the  undisturbed  magnetic 
axis  are  preserved,  and  the  mapping  approximates  the  Identity 
mapping  In  a  neighborhood  of  the  origin. 
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2-3.   The  Form  of  the  Elliptic  Solution 

We  noted  earlier  that  the  fast  wave  envelope  ABCD 

determines  the  boundary  of  the  range  of  influence  of  an 

initial  disturbance  at  (x,y,t)  =  (0,0,0).   For  this  reason, 

we  shall  take  the  solution  (p,u,v,B  ,B  )  to  be  zero  outside 

X  y 

the  fast  wave  envelope.   We  satisfy  the  jump  relations 
across  the  fast  wave  envelope  by  taking  the  solution  to  be 
continuous  across  this  characteristic  surface.   There  are 
three  reasons  for  this  choice.   First,  it  is  the  simplest 
way  to  proceed,  and  if  we  can  construct  such  a  solution,  we 
know  by  our  uniqueness  theorem  that  this  is  the  correct  way 
to  proceed.   Secondly,  Weitzner's  fundamental  solution  [11], 
which  is  a  homogeneous  function  of  (x,y,t)  of  degree  minus 
one,  strongly  suggests  this  choice.   Lastly,  Gardner's 
diffraction  problem  [5]  proceeds  in  this  manner.   Thus,  in 
the  complex  z-plane,  we  take  our  solution  to  be  zero  on  both 
the  upper  and  lower  edges  of  the  slits  from  -A  to  -oo  and 
from  A  to  00  (these  slits  are  the  image  of  the  fast  wave 
envelope  ABCD) . 

We  shall  take  the  solution  to  be  a  constant  state  inside 
the  cusped  portions  of  the  slow  wave  envelope,  and  again 
satisfy  the  jump  relations  by  taking  them  to  be  continuous 
across  these  portions  of  the  slow  wave  envelope.   (The 
constant  state  in  cFaE  need  not  be  identical  with  the  constant 
state  in  GIHK. )   These  restrictions  on  the  form  of  the 
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solution  are  motivated  by  the  same  considerations  as  stated 

In  connection  with  the  nature  of  the  solution  across  the 

fast  wave  envelope.   In  the  complex  z-plane,  we  take  the 

variables  to  be  constants  on  the  silt  from  -a  to  -c  and 

(not  necessarily  the  same)  constants  on  the  silt  from  c  to  a. 

(These  silts  are  the  Image  of  cPaE  and  GIKH.) 

Finally,  there  may  be  a  contact  discontinuity  across 

the  straight  line  segment  of  the  slow  wave  envelope  (Gc) 

which  lies  outside  the  piston.   On  this  line  segment  the  unit 

normal  n  =  (n  ,n  ,n  )  is  given  by 
X  y  t 

n  =  (0,1,0)  (2-3.1) 

Referring  to  the  jump  relations  (1-4.1)  through  (1-4.6)  we 
see  that  we  must  require  that 


a2[p]  +  A^[B^]  =  0  and   [v]  =  [B^l  =  0   (2-3-2) 


across  these  two  line  segments. 

We  now  use  the  representations  (2-1.4)  and  (2-1.5)  in 
the  formal  solution  given  by  equations  (1-3.14)  through 
(1-3.22),  recalling  to  take  the  imaginary  parts  of  these 
Integrals.   We  write  f(T)  for  V(t)  and  have 

z 

P^^^  "  ^^^^"^  \\-~2 — "^2  ~^2 7-   -^f  ('^)d^    (2-3.3) 

J((t  -a  )(t  -c  ) 
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u(z) 


/ac  \  -r 


/     .'   -  a' 


f (T)dT 


(2-3.4) 


v(z]  =  Im 


f (T)dT 


00 


(2-3-5) 


B,{z) 


(^'  I™ 


J 

00 


^-^   2         T 
T    -  C 


(2-3.6) 


B  (z)  =  -Im 


f(Tl 


dr 


J 

00 


(2-3.7) 


The  square  root  in  (2-3.6)  Is  defined  by  equation  (2-1.6) 
The  square  roots  in  (2-3-3)  and  (2-3-4)  are  defined  by  the 
relation 


-^     ;:2" 

T   -  A 


=  (^^ 


"2    2v  ,  2    2T  ~  ^""2    5"^  y    ^^   ? 

T   -a)(T   -Gj     T   -a   I  T-C 


I  2   2   2   2 

(T-a  )(T^-A  )   (2-3.8) 


We  take  f(z)  to  be  analytic  in  the  silt  z-plane, 
figure  9 J  outside  the  image  of  the  piston.   The  conditions  on 
the  slow  and  fast  wave  envelopes  (more  precisely,  on  the 
corresponding  slits)  are  satisfied  by  taking  f(z)  to  be 
real  for  real  z  outside  the  piston  image.   By  oo  in  these 
integrals  we  mean  +00  on  the  real  axis-   To  satisfy  the  condi- 
tions on  the  slits  in  the  left-half  plane  we  must  be  able  to 
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integrate  around  a  large  simlclrcular  path  from  +00  to  -00 
without  picking  up  a  contribution  to  the  various  integrals 
To  Insure  that  we  can  do  this  we  require  that 


f(z) '^  ^  (1  +  0(^))   for  large  |z|  (2-3-9) 


An  examination  of  the  Laurent  series  of  such  a  function 
(about  the  point  00)  shows  that  it  will  be  real  on  the 
imaginary  axis.   Finally,  we  require  that  any  singularities 
which  f(z)  may  possess  lie  within,  or  on  the  surface  of  the 
image  of  the  piston  in  the  z-plane. 
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3.   A  Solution  Involving  only  the  Complex  Characteristics 

In  this  section  we  shall  use  the  formalism  developed 
in  the  last  section  to  construct  a  solution  depending  only 
on  the  complex  characteristics.   (Other  solutions  which  we 
will  discuss  later  employ  both  the  complex  and  the  real 
characteristics.)   Such  a  solution  will  be  possible  only 
for  special  geometries,  and  we  treat  the  shape  of  the  piston 
as  one  of  the  variables  of  the  problem.   Thus,  we  are  led 
to  consider  a  free  boundary  problem. 

The  free  boundary  problem  is  formulated  in  the  second 
subsection.   In  the  last  subsection  we  carry  out  the  solution 
of  this  problem  by  conformal  mapping.   It  turns  out  that  to 
first  order  in  the  expansion  speed  of  the  piston  the  resulting 
(unique)  piston  cross-section  is  circular.   The  explicit 
solution  to  this  problem  is  exhibited.   As  mentioned  earlier, 
the  linearized  problem  is  not  evidently  a  faithful  represen- 
tation of  the  nonlinear  problem.   The  presence  of  a  magnetic 
boundary  layer  is  indicated. 
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3-1 .   Formulation  of  the  Problem 

In  order  to  find  a  solution  Involving  only  the  complex 
characteristics,  we  must  find  an  analytic  function  f(z) 
which  has  the  properties  discussed  at  the  end  of  the  last 
section,  and  which,  in  addition,  satisfies  the  appropriate 
boundary  conditions  on  the  surface  of  the  piston,  namely 
(l-5-ll)  and  (1-5-12).   We  shall  consider  slow  expansions 
of  the  piston,  so  that  the  profile  in  the  (x-,  -p) -plane 
(and  thus  in  the  z-plane)  will  be  small.   From  (2-2.1)  we 
see  that  to  fourth  order  we  can  write 


^  -  ^'Va^)(1  -/7a^)  =  1   (|z|  small)         (>l.l) 
1  -  z  /c 


Using  this  approximation  for  small  |z|  on  the  image  of  the 
piston  (in  the  z-plane),  and  using  (2-3.4)  and  (2-3.5)  we 
can  write 


u(z)  =  C^  -  a^  Re 


z 

-2^^^  dT  (3-1.2) 

a   -  T 


+ 


IP 


z 
f 


v(z)  =  C   +  Im 
^  '  V 


f(^)dT  (5-1.3) 


±t 

where 


,tp 


C„  =  {^)  im 


2     2 

'^    ~  ^  f(T)dT  (3-1.4) 


J 

00 


^2  2s  ,  2     27 

,T   -  a  )(t   -  c  ) 


-  46 


and. 

+t 

r 


f(T)dT  (3-1.5) 


are  constants.   The  point  p  is  the  intercept  of  the  piston 
with  the  positive  (^)-axis,  while  t  is  the  intercept  with 
the  positive  (|-)-axis  (see  figure  9)-   The  various  plus  and 
minus  signs  allow  for  the  different  quadrants  which  may  be 
under  consideration. 

Since  f(z)  is  to  be  real  for  real  z  (to  satisfy  the 
conditions  at  the  slow  and  fast  wave  envelopes)  and  have  a 
second  order  zero  at  infinity  we  have  immediately  that 


C^  =  0  (3-1.6) 


Furthermore,  as  we  mentioned  earlier,  an  examination  of 
the  Laurent  series  about  the  point  infinity  for  the  function 
f(z)  shows  that  f(z)  will  also  be  real  for  z  on  the  imaginary 
axis  (outside  the  image  of  the  piston).   Thus,  we  have  that 


C^  =  0  (3-1.7) 


We  now  introduce  arc  length  s,  along  the  image  of  the  piston 
as  shown  in  figure  10.   On  the  image  of  the  piston  we  write 

z(s)  -  ^(s)  +  iii(s)  (3-1.8) 


-  ^7 


If  we  write 


e(z(s))  =  0^(^(s),  n(s))  +  10.(^(s),ii(s))        (3-1.9) 


we  have  that 
z(s) 


e(T)dT  = 


r 


(^(a),ri(a)) 


+  i0^(4(a),r,(o))j  U^{o)   +   ±']^{o)\   da    (3-1.10) 


where  s   is  the  value  of  arc  length  corresponding  to  the 
point  z  .   Using  (3-1-10)  we  concliide  that 


z(s) 


d 
ds 


Im 


(T)dT';  -  e^(^(s),Tl(3))TlJs) 


+  9^ii{^),n{s))i^{s) 


z(s) 


ds 


Re 


:T)dT^  =  0^(^(5), n(s))e^(s) 


; 


■3-l.il 


-  9^{i{s),n{s))n^{s) 


(3-1.12) 


For  a  point  z(s)  on  the  image  of  the  piston  we  can  now  vrrJte 
(3-1.2)  and  (3-1-3)  in  the  form 

z(s) 


u(z(s)) 


2  ^ 
a  Re 


Illl 


o2     2 
a   -  T 


dT 


(3-1.13) 
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e(s) 
v(z(s))  =  Im  I   f(T)dt  (3-1.14) 

P 

Now,  since  the  z-plane  is  the  image  of  the  (x-,  x-) -plane 
under  the  mapping  (2-1.8)  we  have  that 

dU. 

di^  (P)  =  ?JP)  0-1-15) 

and 

dUo 

^  (P)  =  nJP)  (3-1.16) 


ds   ^  ^    's 

where  P  is  a  point  on  the  image  of  the  piston.   Differentiating 
(3-1.13)  and  (5-1. 1^),  with  the  aid  of  (3-1.11)  and  (3-1.12), 
we  find  that  (3-1.15)  and  (3-1-16)  are  equivalent  to  the  two 
conditions 

-a2p^(|(s),ri(s))|Js)  +  a2F.(^(s),Ti(s))iiJs)  =  ^Js) 

(3-1.17) 

fr(^(s),r)(s))rijs)  +  f.(e(s),Ti(s))^Js)  =iijs)   (3-1.18') 
where  we  have  written 


F(z(s))  =   ^^^^1^^   =  F  (|(s),n(s))  +  iF.(^(s),n(s)) 
a  - z  (s) 


(3-1.19) 
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I     1 2 
To   order     z       we  write 


F(,(3))   .  IM£I1  (5-1.20) 

a 

so  that  (3-1-15)  and  (3-1.16)  become 

-fr(^(s),n(s))e3(s)+f.(e(s),Ti(s))n3(s)==^g(s)    (3-1.21) 
f^(^(s),ri(s))n3(s)  +  f.(^(s),ii(s))^3(s)  =113(3)    (3-1.22) 

In  these  relations,  we  regard  not  only  f(z(s))  but  also 
z(s)  as  unknown  quantities.   Thus,  to  find  the  allowable 
piston  profiles  which  will  yield  flows  depending  only  on 
the  complex  characteristics  we  must  solve  a  free  boundary 
problem.   We  shall  formulate  this  free  boundary  problem  in 
the  next  subsection. 


3-2 .   Formulation  of  the  Free  Boundary  Problem 

We  denote  by  S  the  slope  of  the  yet-to-be  determined 
image  of  the  piston  cross-section.   That  is 

n  (s) 
S(s)  =  ,^j^  (3-2.1) 

s 

We  can  rewrite  equations  (3-1.21)  and  (3-1.22)  in  the  form 
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-f^(^(s),ri(s))    +    S(s)f.(^(s),il(s))    =    1  (3-2.2) 


f^(|(s),Ti(s))    +    (gAy)    fi(^(s),ri(s))    =    1  (3-2.3) 


Equation    (5-2.3)    gives 


f.(e(s),i(s)) 


Using  the  result  (3-2.4)  in  (3-2.2)  we  find  that 


f^(e(s),ii(s))  +  f2(^(s),ii(s))  =  1  (3-2.5) 


We  shall  use  equations  (3-2.4)  and  (3-2.5)  to  solve  our 
problem. 


3-3.   The  Solution  of  the  Problem 

We  consider  the  mapping  properties  of  the  function 
^  =  f(z).   This  function  maps  the  region  exterior  to  the 
piston  image  in  the  z-plane  onto  the  unit  disc  in  the  ^ -plane, 
infinity  being  mapped  into  the  origin  (see  figure  11).   As 
the  image  of  the  piston  is  traversed  in  the  counterclockwise 
sense,  the  unit  circle  of  the  i^-plane  is  traversed  in  the 
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clockwise  sense.   If  we  write 


f(z(s))  =  e^^^^)  (3-5.1) 


then  equation  (3-2.5)  Is  automatically  satisfied.   Equation 
(3-2.4)  gives 

Thus  (see  figure  12)  we  have 

tan  -^^(3)  -  cot(^^)  =  tan(5  -  ^^)  (3-3.3) 

so  that 

^^(s)  =  f  -  ^-^  +  riTT        (n  an  Integer)       (3-3-4) 

From  the  geometry  of  figure  12  It  Is  seen  that  the  appropriate 
choice  Is  n  =  0,  so  that  (3-3.4)  gives 

^(3)=|-^  (3-3.5) 

Referring  to  equation  (3-3.1)  we  see  that 

^(s)  =  arg  f  (5-3.6) 

Equation  (3-3-5)  gives 
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arg  f  =  IT  -  2V(s)  (3-3.7) 

We  will  now  compute  the  change  in  arg  f  as  we  traverse  the 
piston  image.   It  is  convenient  to  start  at  the  bottom  of 
the  piston  image  (see  figure  12)  where  i/^  =  0  =  27r  and  proceed 
counterclockwise.   (The  corresponding  value  of  P/2    is  ^ir/2 , 
and  the  unit  circle  is  traversed  in  the  clockwise  sense.) 
We  find  that 

A  arg  f  =  -  47r  (3-3-8) 

We  have  traversed  the  piston  image  in  the  negative  sense 
(we  are  considering  the  domain  exterior  to  the  piston  image) 
and  so  we  have,  by  the  argument  principle,  that 

N  -  P  -  2  (3-3.9) 

where  N  is  the  number  of  zeros  of  f  in  the  domain  exterior  to 
the  piston  image,  P  being  the  corresponding  number  of  poles. 
Since  the  function  f  is  to  be  regular  in  this  exterior 
domain,  P  =  0,  and  then  (3-3-9)  gives  N  =  2.   Equation  (2-5.9) 
requires  that  f  have  a  double  zero  at  infinity.   Thus,  we  can 
write 


f(z)  =  ^  (3-3.10) 

z 

where  F(z)  is  regular  in  the  exterior  domain  of  the  z-plane 
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and  free  of  zeros.    Now,  we  can  Introduce  the  square  root 
of  the  analytic  function  f(z),  without  needing  a  branch  cut. 
We  write 

g(z)  ^/rr^  ./IpT  (3-3.11) 

Equations  (3-2.4)  and  (3-2.5).  with  the  aid  of  (3-3-1)  and 
(3-3.2)  give 

S^(^(s),n(s)) 
s(^)  =  g>(s),n(3))  ^>5-^2) 

g^(^(s),n(s))  +  g^(^(s),Ti(s))  =  1  (3-3.13) 

The  function  g(z)  has  a  first  order  zero  at  infinity  in 
the  z-plane.   The  point  infinity  in  the  plane,  under  the 
mapping  w  -   g(z),  maps  into  the  origin  in  the  w-plane. 
We  denote  the  inverse  mapping  by 

z  -  k(w)  (3-3-14) 

The  function  k:(w)  will  have  a  first  order  pole  at  the  origin 

(see  figure  13)-   We  now  exploit  one  of  the  principles  of 

conformal  mapping:   if  the  mapping  z  =  Z(w)  maps  a  curve  "P 

w 

in  the  w-plane  onto  a  curve  P   in  the  z-plane,  and  the 

tangent  to  P  makes  an  angle  y_  with  the  real  axis  in  the 

w  w 

w-plane,  while  the  tangent  to  P*  makes  an  angle  y     with  the 

^  z 

real  axis  in  the  z-plane,  then 
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arg[^||^]  =  Y,  -  Y^  (>5.15) 


Using  (5-3-5)  we  have 

arg  k'(w)=(5-^)  -  (^-5)=Tr-i?5(s)   (>5.l6) 


'2    2^    V  2     2 


Now, 


arg  w  =  ^ 


=^  (3-3.17) 


so  that 


2 
arg  w  =  j6{s)  (3-3-18) 


Adding  (3-3-16)  and  (3-3-18)  we  find  that 


2 
arg  w  +  arg  k'(w)  =  ir  (3-3-19) 


Therefore 

2 
arg[w  k'(w)]  =  (2m  +  l)v        (m  an  integer)   (3-3-20) 

2 
We  draw  the  conclusion  that  w  k'(w)  is  negative  real  on  the 

unit  circle  in  the  w-plane. 

As  we  noted  above,  k(w)  will  have  a  first  order  pole  at 

the  origin,  and  be  regular  elsewhere  inside  and  on  the  unit 

circle.   The  derivative  k'(w)  will  have  the  same  regularity 

properties,  but  will  have  a  second  order  pole  at  the  origin. 

2 
Thus,  w  k' (w)  will  be  regular  on  the  unit  disc  in  the  w-plane. 


-  55  - 


2 

and  real  on  the  unit  circle.   We  can  continue  w  k'(w)  by 

the  Schwarz  reflection  principle  to  a  function  that  is 


regular  in  the  entire  w-plane.   By  Liouvllle's  theorem, 

2 
w  k' (w)  must  be  a  constant,  in  fact,  a  nej 

by  virtue  of  (3-3.20).   Thus  we  may  write 


2 
w  k' (w)  must  be  a  constant,  in  fact,  a  negative  real  constant 


w  k'(w)  =  -  r    (r  a  positive  constant)   (3-3-21) 


If  we  solve  this  differential  equation  for  the  function 
k(w)  we  find  that 

k(w)  =  ^  (3-3.22) 


By  equation  (3-3-1^)  we  finally  have 


-  =  ^  05.23) 


Multiplying  equation  (3-3-23)  by  its  complex  conjugate  we 

find  that 

2 
|z|^  =  ^^  (3-3-24) 


w 


The  inverse  mapping  z  =  k(w)  maps  the  unit  circle  in  the 

w-plane  onto  the  image  of  the  piston  in  the  z-plane.   Setting 

2 
lw|   =  1  in  equation  (3-3.24)  we  find  that 

2     ? 
|z|   =  r^  (3-3-25) 
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which  Is  the  equation  of  the  image  of  the  piston  in  the 
z-plane.   From  (3-3.23),  and  the  fact  that  the  z-plane  is 
mapped  into  the  w-plane  by  the  function  w  =  g(z),  we  find 
that 


g(z)  =  -TT 


(3-3.26) 


Finally,  with  the  aid  of  equation  (3-3-11)  we  find  that 


f(z) 


o 


(3-3.27) 


which  obviously  satisfies  all  the  desired  conditions. 

We  thus  have  that  for  a  sufficiently  small  circular 
piston  of  radius  r  the  solution  is  given  by 


cr 


P(z)  =   (^x-)  i"i 


aA 


LZZZIZZ—di 

~r2         2w    2         27    T 
(t    -  a    )(t     -  c    ) 


(3-3.28) 


u(z) 


acr. 


A 


-)    Im 


2       „2 
T      -A 


J 

00 


/    2         2^  ,    2         2x 
(t     -  a    ) (t     -  c    )    t 


dT 


(3-3.29) 


v(z) 


(r^)    Im 


dr 
2 


y 

00 


(3-3.30) 


cr. 


B,(z) 


(5]r'  I- 


^72         2v  ,    2       ^27    , 
(t     -  a    ) (t     -  a    )   dT 

2         2  3 

T      -  C  T 


(3-3.31) 
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By(z)  =  (-r^)  Im 


z 
r 


00 


where 

It  is  Interesting  to  note  that  the  solution  is  "scaled" 
according  to  the  size  of  the  piston,  but  is  otherwise  inde- 
pendent of  the  shape.   Near  the  piston,  the  magnetic  field 
perturbations  are  of  0(l).   This  indicates  a  magnetic  boundary 
layer.   The  field  near  the  surface  of  the  piston  is  strongly 
distorted  so  that  it  will  be  tangent  to  the  surface  of  the 
piston.   This  feature  of  the  solution  is  not  known  a  priori 
from  the  linearized  equations  or  boundary  conditions. 
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4 .   The  Hyperbolic  Wave 

In  this  section  we  discuss  the  hyperbolic  wave.   This 
part  of  the  solution  Is  constructed  using  the  two  real 
characteristics  from  the  piston  to  the  characteristic  conoid. 
The  first  subsection  discusses  the  form  of  the  solution.   This 
subsection  Is  a  straightforward  extension  of  the  corresponding 
section  which  discussed  the  elliptic  solution. 

The  second  subsection  examines  the  nature  of  this  solution. 
In  particular,  we  point  out  why  a  solution  involving  only  the 
real  characteristics  is  not  possible.   The  final  subsection 
discusses  the  dependence  of  /; ,  the  :r-intercept  of  a  real 
characteristic,  upon  the  variables  j-   and  y.      This  result  will 
play  an  important  role  in  some  later  investigations. 


4-1.   The  Form  of  the  Solution 

In  analogy  with  (2-1.4)  and  (2-1.5)  we  write  for  the 
reciprocal  ^  and  ^   Intercepts  of  a  real  characteristic 


A(0  =  i  (4-1.1) 


-(c)  ~-  (^)  If^^f^y^  (^-1-=) 
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Here  ^  Is  a  real  variable.   Geometrically,  ^  Is  the 
(^) -Intercept  of  any  one  of  the  real  characteristics  from 
the  piston  to  the  characteristic  conoid  (see  figure  14). 
The  two  possible  values  of  ^  are  the  two  real  roots  of  the 
quart ic 

MOil)f  =  [1  -  M0(|)]^        (^-1-5) 

We  shall  assume  that  the  piston  is  symmetric  about  both  the 

^  and  ^  axes.   The  symmetry  about  the  ^r-axis  (the  magnetic 
t     t  ^ 

axis)  is  inherent  in  the  problem,  and  an  essential  require- 
ment for  solutions  to  exist.   The  symmetry  about  the  (^)-axis 
is  assumed  for  the  sake  of  convenience . If  this  symmetry  were 
not  present,  the  lower  limits  in  the  integrals  which  express 
the  solution  would  be  different  in  magnitude  for  the  right 
and  left  hand  hyperbolic  waves. 

We  now  introduce  the  notion  of  a  relative  parity  diagram. 
Such  a  figure  consists  of  a  small  picture  of  the  axes,  with 
the  relative  signs  of  the  function  in  question  placed  in  the 
four  quadrants.   We  arbitrarily  assign  the  positive  sign  to 
the  first  quadrant.   Recalling  that  — yy-y  is  the  (^) -intercept 
of  the  real  characteristic  whose  (ir) -intercept  is  1^ ,   we  can 
immediately  write  down  the  relative  parity  diagram  for  the 
function  M-(C);  it  is 

+   + 
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If  we  now  notice  that  the  factor  ^   i^  (^-1-2)  has  the 
relative  parity 


+ 


+ 


we  can  conclude  that  the  square  root  In  (^-1.2)  must  be 
defined  to  have  the  relative  parity 


+ 


+ 

This  definition  of  the  square  root  serves  to  make  M-(C) 
positive  for  the  real  characteristics  In  the  upper  half 
plane,  and  negative  for  the  characteristics  In  the  lower 
half  plane. 

We  now  define  the  (tt) -Intercept  of  the  real  character- 
istic from  the  top  (or  bottom)  of  the  piston  to  the  right 
hand  cusped  portion  of  the  slow  wave  envelope  to  be  C,      (see 
figure  14),  and  choose  z^  =  +  ^^  In  the  formal  representation 
of  a  conical  solution  given  by  equations  (1-3.14)  through 
(1-3.18).   The  plus  sign  corresponds  to  the  right  hand  hyper- 
bolic wave  while  the  minus  sign  corresponds  to  the  left  hand 
hyperbolic  wave.   As  mentioned  earlier,  if  the  piston  were 
not  symmetric  about  the  (|-)-axls,  we  would  have  to  introduce 

■p  T 

^   and  Cp,  fo^  the  two  sides.   We  can  now  write  the  hyperbolic 
wave  in  the  form 
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P(C) 


(— ) 


2   .2 
T  -  A 


'^^  I  V^^  7^7712  12. 


±C 


(t  -a  )(t  -  c  ) 


Tg(T)dT 


) 


K(f^  |) 


u(0  = 


v(C)  = 


^  A  ^ 


^   f  (t  -a  )(t  -  c  ) 


tC, 


g(T)dT(  K(f,  f) 


+  / 


•t'  t 


(11-1.4) 


g(T)dT    K(J,  f)       (4-1.5) 


(4-1.6) 


/ 


B^(C) 


(— ) 
^aA^ 


tc, 


T^   2772  727 
(t  -  a  ) (t  -  a  J 

2    2 

T   -  C 


■) 


g(T)dT 


'  -^(f'  I 


(4-1.7) 


B^(0 


+ 


iC, 


ill 


X  y. 


dx   K(f,  i) 


(4-1.8) 


where 


=  ?(p  f) 


(4-1.9) 


The  plus  signs  In  the  lower  limits  of  integration 
refer  to  the  right  half  plane,  while  the  minus  signs  refer 
to  the  left  half  plane.   The  function  K(^,  |^)  is  defined  to 
he  +1  in  the  shaded  regions  of  figure  15  and  zero  elsewhere. 
Figure  l6  is  an  enlargement  of  this  shaded  region  around  the 
cusp  of  the  slow  wave  envelope  which  is  on  the  positive 
(:r-)-axis.   The  function  g(C)  is  a  real  function.   We  take 
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the  relative  parity  of  this  function  to  be 


+ 


+ 

Heuristically,  we  can  argue  that  this  is  the  correct 
choice.   It  is  reasonable  to  suspect  that  v  will  be  positiv^^ 
in  the  upper  half  plane  and  negative  in  the  lower  half  plane 
An  examination  of  (4-1.6)  shows  that  our  choice  for  the 
relative  parity  of  g(C)  achieves  this,  provided  that  the 
actual  sign  of  g(^)  in  the  first  quadrant  is  negative.   In 
the  next  section  this  choice  of  parity  will  be  rigorously 
justified . 


4-2 .   The  Nature  of  the  Solution 

Each  characteristic  from  the  piston  to  the  cusped  portion 
of  the  slow  wave  envelope  determines  a  value  of  ^.   Thus,  the 
solution  we  have  written  down  is  constant  along  each  of  the 
real  characteristics.   In  several  respects  it  is -analogous 
to  a  simple  wave  solution.   By  choosing  the  lower  limits  of 
Integration  in  (4-1.4)  through  (4-1.8)  to  be  +  Cq  we  have 
insured  that  the  variables  are  zero  on  the  straight  line  portions 
of  the  boundary  of  the  hyperbolic  wave;  consequently,  no  prob- 
lems of  discontinuity  arise  on  these  portions  of  the  boundary 
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of  the  hyperbolic  wave  (see  figure  l6). 

As  mentioned  above,  each  of  the  variables  is  constant 
on  a  given  characteristic.   Actually,  such  a  characteristic 
extends  to  infinity,  but  we  do  not  want  any  of  the  variables 
to  be  nonzero  outside  of  the  fast  wave  envelope,  since  this 
envelope  determines  the  range  of  Influence.   Thus,  at  some 
point  on  each  characteristic  the  values  of  the  variables 
must  abruptly  change  from  constants  to  zero,  i.e.,  there  will 
be  discontinuities  in  the  hyperbolic  wave.   These  discontin- 
uities must  occur  on  a  characteristic  surface.   The  only 
point  inside  the  fast  wave  envelope  where  the  real  character- 
istics are  tangent  to  a  characteristic  surface  is  where  they 
meet  the  slow  wave  envelope.   For  this  reason  we  have 
introduced  the  factor  K(ir,  y)   which  makes  all  the  variables 
zero  beyond  the  point  of  tangency  on  the  cusped  portion  of 
the  slow  wave  envelope. 

Across  that  portion  of  the  slow  wave  envelope  which  is 
shown  in  figure  l6  (but  not  dotted)  the  hyperbolic  wave  is 
discontinuous,  and  does  not  satisfy  the  jump  relations.   It 
Is  for  this  reason  that  no  solution  involving  only  the  real 
characteristics  is  possible. 


-  64 


^-3-   The  Dependence  of  ^  on  ^  and  ^ 


We  know  from  the  first  section  that  A  and  |i.  satisfy 
the  following  two  equations: 

t  -  A(c)x  -  ^L(C)y  -  0  (^-3.1) 


1  -  (a^  +  A^)(a2(0  +  li^iO) 

+  a.^A^-h^{ai->^^{a   +  ^L^(C))  =  0      (4-3.2) 


Equation  (4-3.1)  gives  A((^)  as  a  function  of  M-(C)  and  x- 
and  x-j  namely 


MO  4=^  -^(^'^  (*-3.3) 


where  we  have  written 


X  =  ^  (4-3.4) 


Y  =  f  (4-3.5) 


Inserting  the  representation  (4-3.3)  Into  equation  (4-3.2), 
we  obtain,  after  some  manipulations,  an  equation  relating  i^ 
and  Y   and  ^  : 

(Y^-  (a^  +  A^))C^+2(a^  +  A^)XC^  +  (a^A^  -  (a^  +  A^  )  (X^  +  Y^  )  )C^ 
-  2a^A^XC  +  a^A^(X^  +  y^ )  ^  0  (4-3-6) 
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where 

C  =   C(x,Y)  (^-3.7) 

We  now  obtain  i^  as  a  power  series  In  X  and  Y  by 
employing  the  two-dimensional  form  of  Taylor's  theorem, 
namely 


C(X,Y)  =  c(o,o)  +  Cx(0'0)X  +  ^y(o,o)y  (4-3-8) 


By  implicit  differentiation  of  equation  (4-3.6)  we  find 
that 

IJ  (0,0)  =  ±  c  (4-3.9) 

II  (0,0)  =.  0  (4-3.10) 

^   (0,0)  =:  0  (4-3.11) 

hX 

^2 

^xSy  ^°'°^  "  °  (4-3.12) 

n2  3 

f4  (0,0)  =  +  -^  (4-3.13) 

Sy^  a'^A 

Thus,  -using  equations  (4-3.9)  through  (4-3.13)  in  equation 
(4-3.8)  we  obtain 
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C(X,Y)  =  +  c[l  +  ^  p  p  +  .. .] 

2  a'' A 


(^-5.1^) 


or,  using  (4-5.4)  and  (4-3-5),  we  have 


C(J,  f)  =t  C[l4--f^  (f)   + 


•t'  t 


]        (4-3.15) 


2a  A 
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5-   The  Solution  for  a  Symmetric  Piston  with  a  Full 
Hyperbolic  Wave 

In  this  section  we  shall  construct  the  solution  for  an 
arbitrary  piston  whose  cross-section  is  symmetric  about 
both  the  (^)  and  (|^)-axes,  although  we  briefly  indicate  how 
the  solution  might  be  obtained  for  a  piston  sjnnmetric  about 
the  (^)-axis  only,  when  such  a  solution  is  possible.   We 
further  suppose  that  the  entire  surface  is  covered  by  the 
hyperbolic  wave;  this  simply  says  that  the  top  and  bottom 
of  the  piston  are  not  too  flat.   This  latter  restriction  will 
be  dropped  in  the  next  section.   The  solution  is  constructed 
by  employing  both  the  real  and  complex  characteristics. 

In  the  first  subsection  we  write  down  the  formal  solution 
which  we  shall  use.   The  second  subsection  discusses  the 
boundary  conditions  which  we  shall  impose  on  this  formal 
solution.   Roughly  speaking,  there  are  two  conditions  which 
are  not  Immediately  satisfied;  these  are  the  conditions  on 
part  of  the  slow  wave  envelope,  and  the  conditions  on  the 
piston.   The  third  and  fourth  subsections  discuss  how  we 
satisfy  these  remaining  conditions. 

The  next  three  subsections  discuss  the  solution  of  the 
relevant  boundary  value  problem.   Various  singularities  may 
appear  in  the  solution,  and  the  nature  of  these  singularities 
is  examined.   Subsection  eight  discusses  the  determination 
of  the  one  arbitrary  constant  which  occurs  in  the  solution. 
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The  section  ends  with  a  discussion  of  the  possible  extension 
of  our  methods  to  a  piston  whose  cross-section  Is  not  sym- 
metric about  the  (^)-axls. 


5-1.   The  Form  of  the  Solution 


We  write  our  solution  as  the  sum  of  an  elliptic  solution 
and  a  hyperbolic  wave.   Using  equations  (2-3.3)  through 
(2-3.7)  and  equations  (4-1.4)  through  (4-1.8)  we  can  write 
our  formal  solution  in  the  form 

z 
m  \l — - — '^^  "  ^^ —   Tf(T)dT     (5-1.1) 


P  ih    f )  =  (^)  I 


■t'  t 


aA' 


X 


-<(ij) 


tc, 


(t  -  a  )(t  -  c  ) 


^^._J)',^^_^.)  Tg(t)d.JK{f.  £) 


u(f.  |)  -    (^)    in, 


2   „2 
T  -  A 


J 
00 


— ^ 2 2 2"  ^i'^)^^         (5-1.2) 

(t  -a  )(t  -  c  ) 


2   „2 
T  -  A 


(T  -  a  j  (^T  -  c  ) 
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V 


(|,  f^)  =  Im  [  f(T)dT  +    J  g(T)dT  K(|,  f) 

—  ^r\ 


(5-1.3) 


B.(^.  ^)M^)lJp^Iapi?IZliMa.   (5-1.4) 


J 
00 


T   -  C 


+  < 


c  

J   /         T   -  C  I 


By(|.  |)  =  -I- 


£M 


r  / 


dT 


111 


dT 


>.  + 


^L 


K(|,  f)   (5-1.5) 


where 


and 


C  - 


zf-  ^ 


^(^  ^) 


(5-1.6) 


(5-1.7) 


As  before,  f(z)  is  an  analytic  function  of  the  complex 
variable  z,  and  g(C)  is  a  real  function  of  the  real  variable 
C-   The  function  K(:r-,  ^)  Is  again  defined  to  be  +1  In  the 
shaded  region  of  figure  15,  and  zero  elsewhere.   In  the  real 
Integrals,  the  plus  signs  In  the  lower  limits  of  Integration 
refer  to  the  right  half  plane;  the  munus  signs  refer  to  the 
left  half  plane. 

To  satisfy  the  boundary  and  matching  conditions,  which 
we  shall  discuss  In  the  next  subsection,  we  will  place 
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certain  restrictions  on  f  and  the  path  of  Integration  in 
the  contour  integrals.   The  boundary  conditions  on  he 
surface  of  the  piston  serve  to  complete  the  determination 
of  f ( z )  and  g ( C ) • 


5-2 .   The  Boundary  Conditions 

As  before,  we  take  the  solution  to  be  zero  outside  the 
fast  wave  envelope  ABCD,  thus  preserving  the  range  of 
influence.   Proceeding  as  before,  we  take  the  solution  to 
be  continuous  across  that  characteristic  surface.   Because 
of  the  factor  K(^,  y)   which  multiplies  the  real  integrals, 
the  solution  in  the  neighborhood  of  the  fast  wave  envelope, 
ABCD,  consists  only  of  the  contour  Integrals.   We  satisfy 
these  requirements  by  taking  the  analytic  function  f(z)  to 
be  real  on  the  slits  from  -A  to  -oo  and  from  A  to  oo  (see 
figure  9).   The  function  f(z)  must  also  have  a  double  zero 
at  infinity. 

Inside  that  portion  of  the  cusped  slow  wave  envelopes 
which  is  not  covered  by  the  hyperbolic  wave  (see  figure  l6) 
the  solution  still  consists  of  only  the  elliptic  solution. 
Again,  we  require  the  solution  to  be  a  constant  state  in 
these  regions,  and  take  it  to  be  continuous  across  the 
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bounding  portion  of  the  slow  wave  envelope.   This  condition 

is  satisfied  by  taking  f(z)  to  be  real  on  the  real  axis 

from  -A  to  -t      and  from  C   to  A.   Coupled  with  the  fast  wave 
o  o 

envelope  conditions,  we  now  have  that  f(z)  must  be  real  for 
2  real  and  -oo  <  z  <  -^   or  T  <  z  <  oo.   Furthermore 


f(z)  =  ^  (1  +  0(i))    (|z|  large)  (5-2.1) 

z  '  ' 


That  portion  of  the  cusped  slow  wave  envelopes  across 

which  the  hyperbolic  wave  exhibits  a  discontinuity  maps 

(under  (2-1.8))  onto  the  two  silts  from  -T   to  -c  and  from 

^   to  c.   (These  are  really  pieces  of  the  silts  from  -a  to 

-c  and  c  to  a;  see  figure  l?-)   We  shall  take  the  solution 

to  be  continuous  across  this  portion  of  the  slow  envelope 

also,  but  here  both  the  elliptic  and  hyperbolic  portions  of 

the  solution  must  be  considered.   Thus,  we  require  that  the 

elliptic  solution  on  the  slits  from  -f   to  -c  and  from  c  to 

^o 

r   match  the  sum  of  the  hyperbolic  wave  and  the  constant 
state  in  the  other  portions  of  the  slow  wave  envelope.   We 
shall  now  state  this  mathematically. 
The  constant  state  is  given  by 


+ 


C 


(t  -a  )(t  -c  ) 


+  00 
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±^ 


+ 


o 


±00 


(t    -a    )(t   -c    ) 


^c^l'    t)    =   ^"^ 


+  00 


2~2 

2    ^"2         2      ^      f  ('^)dT 
(t    -a    )(t    -c    ) 


(5-2.3) 


(5-2.4) 


+ 


±? 


B;c(|'f)  =  (^)  ^- 


o 


J 

+  00 


iZzSIZzZi  fill  dT 


(5-2.5) 


V^t'  t) 


-    Im 


+00 


fill 

T 


dT 


(5-2.6) 


In  these  formulas  the  plus  signs  refer  to  the  right  half  plane; 
the  minus  signs  refer  to  the  left  half  plane.   It  should  be 
noted  that  the  only  contribution  to  these  integrals  occurs 
from  +A  to  +a.   The  total  solution  will  be  continuous  across 
the  slow  wave  envelope  if  we  satisfy  the  following  five 
conditions: 


Po(f '  f)  +  (ij) 


±c 


/  r'-A' 

,2      2s  f    2      2.    ^ 
(t  -a  )(t  -c  ) 


g(T)d 


T 


+00 


2—^7  '^f  ('^)d'r 

T  -C   ) 


(5-2.7) 
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J 


'  2_   2 

(t   -a    )(t   -c    ) 


o 


(•ac 


(^)    im 


2    .2 
T    -A 


y   ,    2      2w    2      2x 
J      /   (t    -a    )(t    -c    ) 

+00 


f(T)dT  (5-2.8) 


v;(^,   I)   +    \      g('^)dT  =   Im 


f (T)dT 


+  c 


(5-2.9) 


B;e(|'    f)    -    (^) 


/^       /T~2       2x/^2     „2 


T-^-a")(t"-A")    g(T) 


2      2 

T     -C 


dT 


±c, 


==  (^)  ^- 


'^      ^^T^-a^)(T^-A^)    fJV 


2      2 

T     -C 


dT  (5-2.10) 


+  CO 
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R-  r^  y^  - 
V^t'  t^     j 


ill 


dT   =    -   Im 


fill 


dT  (5-2.11) 


+c 


o 


where 


c  =  a|,  ?) 


(5-2.12) 


and 


c  <   kl  <  C 


(5-2.13) 


If  we   now  use   equations    (5-2.2)    through    (5-2.6)   we    see 
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that  the  conditions  (5-2-7)  through  (5-2. 11)  can  be  written 
In  the  form 

c 

&    I  /  2    ^2 "^'^2-?-  '^(s(t)  -  Im  f(T))dT  =  0    (5-2.14) 

^^   J  /  (T^-a^)(T^-C^) 

^o 


(— ) 

^  A  ^ 


2  a2 
T  -A 


+ 


/  2   2v  /  2   2s 
(t  -a  ) (t  -c  ) 


(g(T)  -  Im  f(T))dT  =  0    (5-2.15) 


n, 


±C, 


■  c  ■ 
•aA' 


:(t)  -  Im  f(T))dT  =  0 


C  

(       /7"2   2T7~2~2T 

/       T  -C 


(5-2.16) 


Im  f(T))dT  =  0   (5-2.17) 


*C, 


J 


(g(T)  -  Im  f(T))  ^  =  0 


(5-2.18) 


We  can  satisfy  equations  (5-2.14)  through  (5-2.18)  by 
taking 


Xi)    =  Im  f{i)  c  <  1^1  <  C, 


(5-2.19) 


We  shall  discuss  this  equation  in  some  detail  in  the  next 
Subsection. 

On  that  portion  of  the  line  segment  Gc  which  lies 
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outside  the  piston  we  again  allow  for  a  contact  discon- 
tinuity, by  requiring  that 


a^[p]  +  A^[B^]  =0    [v]  -  [B^  1  -  0         (5-2.20) 


across  those  line  segments.   In  connection  with  this  possible 
contact  discontinuity,  we  restrict  the  path  of  Integration 
In  the  contour  Integrals  to  cross  the  real  axis  only  In  the 
two  Intervals  (-A,-a)  and  (a, A). 

Finally,  we  require  that  at  a  point  P  on  the  surface 
of  the  piston 

u(P)  =  U^(P)  (5-2.21) 

v(P)  =  U2(P)  (5-2.22) 

To  sum  up,  we  require  that  f(z)  be  real  for  z   real  and 
Cq  ^  1^1  <  °°>  3-"d  have  a  double  zero  at  Infinity.   The  path 
of  Integration  In  the  contour  Integrals  Is  restricted  as 
mentioned  above.   Finally,  we  require  that  (5-2.19)  through 
(5-2.22)  be  satisfied. 
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5-3-   The  Matching  on  the  Slow  Wave  Envelope 

In  this  subsection  we  will  discuss  how  we  satisfy 
equation  (5-2.19).   We  will  denote  the  order  of  the  expan- 
sion speed  of  the  piston  by  e;  more  precisely,  the  maximum 
of  the  y-component  of  the  expansion  velocity  of  the  piston 
is  denoted  by  e.   Then,  by  (4-3.15)  we  have 


^o  =  c  +  0(6^)  (5-3.1) 


2 
Thus,  the  length  of  the  interval  (c,Cq)  is  of  order  e  .   The 

solution  to  our  problem  is  nonuniform  in  e.      In  other  words, 

the  problem  is  inherently  of  lower  order  in  the  neighborhood 

of  the  tips  of  the  cusped  portions  of  the  slow  wave  envelope. 

Elsewhere,  the  first  order  problem  is  of  the  same  character 

as  the  second  oraer  problem,  but  in  these  small  regions 

the  first  order  problem  is  of  an  entirely  different  nature 

than  the  set^ond  ord-er  problem.   Thus,  in  these  small  regions, 

we  shall  carry  the  solution  to  second  order  in  the  piston 

velocity.   Recalling  the  relative  parity  of  the  function  g{ll,) , 

we  can  see  from  equation  (5-2.19)  that  it  will  be  helpful  to 

v;rite 


r   ° 


f(-)  hI 


V 


^^d^\F{z)  (5-3.2) 


The  expression  inside  the  bracket  can  be  written  as  the 
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sum  of  two  singular  Integrals.   If  we  take  r(^)  to  be  an 
even  real  function,  we  have 


c 


o 


2 

■JT 


1 
IT 


r?o 

f^d? 

1 

^-Z 

TT 

-c 


r(l) 


Ai 


(5-3.3) 


-C 


Thus,  the  factor  in  the  brackets  will  take  care  of  the  dis- 
continuities of  f(z)  across  the  slits  from  -t,      to  c  and 
from  c  to  (^  .   Furthermore,  for  large  |z|,  this  factor 
behaves  like  l/z  ,  so  we  need  only  require  that  the  analytic 
function  F(z)  be  bounded  at  infinity.   Thus,  in  (5-3-2), 
P(z)  is  an  analytic  function,  outside  the  piston  image,  which 
is  bounded  at  infinity  and  real  for  real  z.   The  function  r 
will  be  chosen  to  satisfy  equation  (5-2.19)-   For  future 
convenience,  we  normalize  the  function  r  by  setting 

-^o 

f  (5-3.4) 


rli) 


We  now  examine  these  singular  integrals,  making  use  of 
the  Plemelj  formulas.   We  denote  the  upper  half  plane  limit 
of  a  function  by  a  superscript  plus,  appending  the  super- 
script minus  to  denote  the  lower  half  plane  limit.   Let 


R, (z)  =  1 

1  ^   ^     IT 


'O 


%  -  z 


(5-3.5) 


For  z   =   1^^,   with  c  <  C^  <  ^^,    the  Plemelj  formulas  give 
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Rt(q)  =  ir(c,)  +  f  j    fif^de 


(5-3.6) 


and 


P 


Ri(Cl)    =    -ir(Ci)   +? 


rlil 


d^ 


(5-3.7) 


where  the  P  preceding  an  Integral  Indicates  that  the  integral 
is  to  be  interpreted  in  the  Cauchy  principal  value  sense. 
Similarly,  if  we  define 

,-c 
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2,  -  z 


-c, 


we  have  for  z  =  Co  ?  "C^  '^  ^o  '^    ~^ '    that 
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^2VS2 
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-C 
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(5-3.8) 


(5-3.9) 


(5-3.10) 


Thus,  bearing  in  mind  the  relative  parity  of  the 
function  g(C)j  S-^*^  noticing  that  the  principal  value  integrals 
are  real,  we  can  satisfy  equation  (5-2.19)  by  requiring  that 


,+ 


r(C)F(C)  =  g(C  )    c  <  Id  <  C 


o 


(5-3.11) 


If  F  and  g  are  determined  (as  they  will  be  by  the  conditions 
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(5-2.21)  and  (5-2.22)  on  the  surface  of  the  piston)  equation 
(5-5.11)  determines  the  function  r. 


5-4.   The  Boundary  Conditions  on  the  Piston 


In  this  subsection  we  discuss  the  conditions  (5-2.21) 
and  (5-2.22).   As  In  subsection  (3-1)  we  write 
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f  (T)dT  + 
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g(T)dT 


(5-4.2) 


where  s  Is  arc  length  along  the  Image  of  the  piston  In  the 

z -plane  (see  figure  10)  and  C   and  C   are  constants.   We 
^       °       '       u      V 

also  have 


f(z) 


v 


i   -z 


(5-4.3) 


where  the  function  r  has  been  normalized  so  that 
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2 
IT 


/o 


v{i) 


d^  =  1 


(5-^.^) 


In  the  complex  Integrals  along  the  image  of  the  piston  we 
have  that 


z^  =  O(e^) 


(5-^.5) 


so  that  in  these  integrals  we  can  take 


f(z)  = 


TT 


r(^ 


di  \   F(z)  =  F(z) 


(5-^.6) 


Thus  we  can  rewrite  equations  (5-4.1)  and  (5-4-. 2)  in  terms 
of  the  analytic  function  F(z),  as  follows: 

z(s) 


^ac 


u(s)  ^  C^  4-  (^)  Im 


tt 


(,^ inke^)  "'^"^ 


t(s) 


.ac 


.  {^) 


+ 


2   2,,  2   2,  S(T)dT 
,T  -a  )(t  -c  ) 


(5-4.7) 


tc, 


v(s)  =  C„  +  Im 


z(s) 


C(s) 


F(T)dT  + 


g(T:)dT 


(5-4.8) 


ip 


tc, 


We  now  differentiate,  and  employ  equations  (5-2.21) 
and  (5-2.22)  so  that  we  have 


S(^)=^s(^) 


(5-4.9) 
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^(s)=^^is) 


(5-4.10) 


or,  using  the  representations  (5-4-7)  and  (5-4.8)  we  can 
write 


^3(3)  = 


(^)  Re) 


z  -A 


2   ^w  2   2.  ^(^)\  %^< 


-a  )(z  -c  ) 


.  c-t)  M  y 


./     2T2 
z  -A 


^2   2,,  2   2v 
(z  -a  )(z  -c  ) 


F(z) 


^3(3) 


+ 


/     2   2 

C  -A 


^^^  y(^2_^2)(^2_^2.  ^s 


C.(s)g(c) 


(5-4.11) 


Ti^(s)  =  Re|F(z)'pg(s)  +  Im(F(z)j^g(s)  +  eg(s)g(C)  (5-4.12 


where 


z  =  z (s) 


(5-4.13) 


and 


C  =  C(3) 


(5-4.14) 


We  shall  now  eliminate  g(C(3))  from  equations (5-4 . 11 ) 
and  (5-4.12)  to  obtain  a  linear  relationship  connecting 
the  boundary  values  of  the  real  and  imaginary  part  of  F(z). 
We  first  introduce  some  notation: 


F(z(s))  =  P^(Us),n(s))  +  iF^(e(s),ri(s: 


(5-4.15) 
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6(x)^  (x)/r^'2:fg  2,  (5-*. 16) 

/  (t    -a    )(t    -c    ) 
5(z(s))    =    6^(^(s),ri(s))    +   15.(^(s),ri(s))  (5-^.17) 


A(t)   =    (^)  'Ji-'^'-^nf-'    I  (5-4.18) 


^ac^    /  7:7- 


Equations    (5-4. 11)    and    (5-4.12)    can  then  be  written   In  the 
form 


^s   =    ^Ms  +   Sl^s)^r   -    (^l^s    -  Ms)Pl  +  ^"'Csg      (5-4.19) 


TI3    =    Ti^F^  +    e^F.    +    ^^g  (5-4.20) 


Finally,  multiplying  equation  (5-4.19)  by  A  and  subtracting 
equation  (5-^.20)  we  obtain  the  desired  linear  equation 
relating  the  boundary  values  F  (^  (s  )  ,r]  (s  ) )  and  F^  (4  (s  )  ,r)  (s  )  ) : 

[  (A5^-l)ng  +  A5.^  JF^  +  [  (A6^-l)e3  -  A5.Tig]P.  =  A^^  -  ^^ 

(5-4.21) 

As  we  did  In  the  third  section,  we  use  the  second-order 
appr oxlmat 1 on 

6(z(s))  =  -1  (5-4.22) 

This  enables  us  to  write  equation  (5-4.21)  In  the  form 
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ct(s)F  (s)  +  p(s)F^(s)  =  y(s)  (5-^.25) 


where 


a 


(s)  =  A(C(s))£  (s)  +  n  fs)  (5-^.24) 


s 


p(s)  =  ^g(s)  -  A(C(s))iiJs)  (5-4.25) 

7(3)  =  n^(s)  -  A(C(s))^Js)  (5-4.26) 


Before  discussing  the  solution  of  (5-4.23)  we  devote  a  sub- 
section to  the  discussion  of  the  relative  parity  of  the 
terms  in  (5-4.21).   In  the  course  of  this  discussion  we 
verify  our  prior  claim  regarding  the  relative  parity  of 
the  function  g(C)- 


5-5  The  Relative  Parity  of  Equation  (5-4.21) 

In  this  subsection  we  will  show  that  the  relative  parity 
of  equation  (5-4.21)  is  such  that  we  can  solve  our  boundary 
value  problem  for  F(z)  in  a  quarter  plane,  and  then  complete 
the  solution  by  two  applications  of  the  Schwarz  reflection 
principle.   In  particular,  this  enables  us  to  construct  F(z) 
so  that  it  will  be  real  for  z  real  or  pure  imaginary,  so 
long  as  z  lies  outside  of  the  image  of  the  piston. 


-  84  - 


In  section  four  we  found  the  relative  parities  of 
A(^)  and  M-(0  could  be  expressed  by  the  following  parity 
diagrams: 


MO 


+ 


^^(C) 


+ 


+ 


+ 


Thus,  since 


A(0  =  ^  [1  -  a2(a2(o  +  ^i2(c))] 


-1 


(5-5.1) 


we  conclude  that  the  relative  parity  of  A(0  is  given  by 
the  diagram 


A(0: 


+ 


+ 


Recalling  our  definition  of  arc  length  (see  figure  10) 
we  have  immediately  that  the  relative  parities  of  £  (s)  and 
ri^(s)  are 


i,{B): 


+ 


+ 


%{s): 

+ 

+ 

+ 

+ 

Thus,  the  relative  parity  of 


^(s)  =.  A(a3))^^(s)  -  11  (s) 


(5-5.2) 


can  be  represented  by  the  diagram 
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Y(s): 


+ 


+ 


+ 


+ 


A  straightforward  examination  of 


^^^)  =  ^  ^^  -  A^(^^^^)  +  ^^(^))^ 


(5-5.3) 


discloses  that  the  real  and  Imaginary  parts  of  6  have  the 
relative  parities 


+ 


+ 


+ 


+ 


+ 


+ 


Thus,  the  expressions 


and 


a  =  (A6   -  l)r|   +  A5.^ 
^   r      s     IS 


g  =  (A5^  -  1)^^  -  A6.7i^ 


;5-5.^) 


(5-5-5) 


have  the  relative   parities  expressed  by  the  following  parity 
diagrams : 


a(s) 


+ 


+ 


+ 


P(s) 


+ 


+ 


We  can  now  write  equation  (5-4.21)  in  the  form 


S(s)F^(s)  +  p(s)F^(s)  =  y(s) 


(5-5.6) 
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From  this  equation,  we  can  deduce  that  the  boundary  values 
of  the  real  and  imaginary  parts  of  the  analytic  function 
F(z)  have  the  relative  parities 


Fp(s): 


+ 


+ 


+ 


+ 


P.(s) 


+ 


+ 


Finally,  referring  to  equation  (5-^.20),  and  noting 
that  1^  (s)  has  the  relative  parity 


Cjs): 


+ 


we  can  infer  that  g(C)  has  the  relative  parity  stated 
earlier,  i.e.. 


:(0 


+ 


+ 


Thus,  if  we  construct  an  analytic  function  P(z)  whose 
boundary  values  satisfy  equation  (5—^.21)  or  (5— 4-. 23)  in 
the  first  quadrant,  and  extend  this  function  by  first 
reflecting  across  the  positive  imaginary  axis,  and  then 
across  the  entire  real  axis,  the  resulting  boundary  values 
will  satisfy  (5-''+.2l)  or  (5-'4-.23)  on  the  whole  image  of  the 
piston.   This  is  the  strategy  we  will  employ. 
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5-6.   Solution  of  the  Boundary  Value  Problem;  Formal 
Discussion 

In  this  subsection  we  discuss  the  construction  of  the 
analytic  function  F(z).   As  suggested  in  the  last  subsec- 
tion we  shall  solve  a  quarter  plane  problem.   On  the  piston, 
we  satisfy  equation  (5— 4-. 23).   On  the  real  and  imaginary 
axes  we  require  F(z)  to  be  real.   This  will  make  C   in 
(5—^-7)  zero,  and  0   in  (5-^-8)  of  order  e  •   It  will  also 
permit  us  to  extend  the  analytic  function  F(z)  to  the  whole 
plane,   by  applying  the  reflection  principle  twice. 

We  begin  by  multiplying  equation  (5— ^-23)  by  a  yet-to-be 
determined  real  function  H(s): 

H(s)a(s)F^(s)  +  H(s)p(s)F^(s)  =  H(s)y(s)    (5-6.1) 

We  regard  the  left  hand  side  of  equation  (5-6. l)  as  the 
boundary  values  of  the  imaginary  part  of  the  product  of 
two  analytic  functions,  the  analytic  function  F(z)  and  some 
other  analytic  function  G(z).   From  (5-6. l)  we  deduce  that 
the  boundary  values  of  the  real  part  of  G(z)  are  H(s)3(s), 
while  those  of  the  imaginary  part  of  G(z)  are  H(s)a(s).   If 
we  write 

h(z)  =  log  G(z)  (5-6.2) 

we  have  that 


Thu 
are 


Im  fh(z)'l  -  arg[G(z)]  (5-6.3) 

s,  the  boundary  values  of  Im  j^h(z)  j  on  the  piston  Image 


Ms)  =  tan-l  (^)  (5-6.4) 


On  the  real  and  imaginary  axes  we  take  the  boundary  values 
of  Im(h(z)^  to  be  zero.   We  thus  consider  the  quarter  plane 
problem  shown  in  figure  l8. 

We  impose  the  additional  restriction  of  boundedness  at 
infinity  to  insure  uniqueness.   In  the  next  subsection  we 
rigorously  discuss  the  solution  of  this  exterior  Dirichlet 
problem,  including  a  discussion  of  the  singularities  which 
may  enter  the  solution.   We  complete  the  determination  of 
the  analytic  function  h(z)  by  employing  the  Cauchy-Riemann 
equations.   The  function   h(z)  will  contain  an  arbitrary  real 
additive  constant.   We  shall  not  carry  this  constant  in  our 
solution  since  it  drops  out  of  the  final  determination  of 
P(z). 

We  can  now  see  that 

h  (s) 
H(s)  =  e  ^  (5-6.5) 

where  h  (s)  denotes  the  boundary  values  of  the  real  part  of 
h(z).   We  consider  the  function 
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P(z)  =  F(z)G(z)  (5-6.6) 

Equation  (5-6.1)  Implies  that  the  boundary  values  of 
Im  |f(2)>  are  given  by 

h  (s) 
k(s)  =  e  "^   y(s)  (5-6.7) 

Again  we  take  the  boundary  values  of  Im{^P(z)>  to  be  zero 
on  both  the  real  and  imaginary  axes,  and  impose  the  condition 
of  boundedness  at  infinity  (see  figure  19)-   We  postpone  the 
discussion  of  possible  singularities  in  the  solution  of  this 
exterior  Dirichlet  problem  until  the  next  subsection. 

We  construct  Re  I  F(z)  >  by  employing  the  Cauchy-Riemann 
equations,  noting  that  this  construction  introduces  an 
arbitrary  real  additive  constant  into  the  function  F(z). 
We  will  write 

F(z)  =  F^(^,Ti)  +  p  +  iF.(^,ri)  (5-6.8) 

where  [""  is  the  arbitrary  constant.   We  will  determine  the 
value  of  this  constant  by  satisfying  the  normalization 
condition  (5-5. '4-). 

We  now  have  determined  the  analytic  function  F(z)  in 
the  first  quadrant;  it  is  given  by 

F(z)  =  F(z)e-^(^)  (5-6.9) 
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By  construction,  this  function  satisfies  (3-6.1)  on  the 
Image  of  the  piston,  and  Is  real  on  both  the  real  and 
Imaginary  axes.   We  extend  this  function  to  the  whole  z-plane 
(exclusive  of  the  Image  of  the  piston)  by  reflecting  across 
the  positive  r|-axls,  and  then  reflection  across  the  entire 
^-axls.   There  may  be  singularities  (introduced  through  dis- 
continuities In  the  boundary  data)  at  the  points  +p  ,  +t 
(see  figure  9).   In  the  next  subsection  we  discuss  these  and 
other  delicate  features  of  the  above  analysis. 


5-7,   Solution  of  the  Boundary  Value  Problem:  Rigorous 
Discussion 

In  this  subsection  we  shall  verify  that  the  procedure 
outlined  In  the  last  subsection  can  actually  be  carried  out, 
and  examine  the  possible  singularities  In  the  solution.   In 
the  boundary  value  problem  for  Im  ^h(z)j  we  Introduced  the 
function 

A(s)  =  tan-1  (||||)  (5-7.1) 

From  (5-^.2^)  and  (5-^.25)  we  have  that 


a(s)  _  A+S  (5-7.2) 

p(¥)  "  1-AS  ^-^    ' 
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where  S  is  the  slope  of  the  image  of  the  piston  in  the 
z-plane. 

By  differentiation  we  find  that  the  maxim-um  value  of  A 
(in  magnitude)  occurs  for 

W-'^'ll  -y^   -2=7a"]  (5-7.3) 

Using  (5-4.18)  it  is  an  easy  matter  to  obtain  the  bound 


A|  <  2{A^    -   c^)  (5-7.4) 


so  that  A  is  bounded  from  above  and  below.   Now  since  S  and  A 
are  nonpositive  in  the  first  quadrant,  being  zero  only  where 
the  image  intersects  the  ^-axis,  we  see  that  the  numerator 
of  (5-7'2)  possesses  only  one  zero.   Thus,  the  data  on  the 
piston  image  all  lie  on  a  single  branch  of  the  arctangent 
function.   We  choose  the  principal  branch 


I  <  tan"^0  <  I  (5-7.5) 


With  this  choice,  the  zero  boundary  values  specified  on  the 
coordinate  axes  are  also  on  the  same  branch  of  the  arctangent 
There  may  be  a  singularity  in  P(z)  at  the  "front"  of  the 
piston,  where  the  boundary  values  can  be  discontinuous. 
Referring  to  the  sketch  of  the  relevant  geometry  (figure  20), 
and  setting  A  =  0  in  this  neighborhood  of  the  ^-axis,  we  find 
from  (5-7-1)  that 
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A^  =.  tan"^S^  (5-7.6) 


o        o 

So  Ions  as 


-cx)  <  Sq  <  0  (5-7.7) 


we  have  that 

-  f  <  A^  <  0  (5-7.8) 

Thus,  for  a  piston  with  zero-angle  cusps  along  the  undisturbed 
magnetic  axis,  the  boundary  values  of  the  harmonic  function 
Im^h(z)'^  will  be  continuous. 

To  discuss  the  general  case,  we  introduce  the  harmonic 
function 

M(S^)arg(z  -  ^^)  (5-7.9) 

where 

tan"  S 

M(S^)  =  if-  (5-7.10) 

IT  +   tan  Sq 

Using  (5-7.8)  we  see  that 

-1  <  m(Sq)  >  0  (5-7.11) 

The  values  of  the  harmonic  function  (5-7.9)  on  the  real  axis 
are  0,  while  on  the  line  segment  of  slope  S  ,  the  values  are 
IT  +   tan"  S  .   Thus,  the  harmonic  function 

lm{h(z)^  -  M(SQ)arg(z  -  ^^)  (5-7.12) 
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assumes  zero  boundary  values  in  the  neighborhood  of  the  "front" 
of  the  piston.   Thus,  locally, 

lm[h(z)]  ~  M(SQ)arg(z  -  ^^)  (5-7.13) 


and 


M(S^) 


h(z)^  log[(z  -  ^^)    "  ]  (5-7.1^) 


Prom  (5-7.14)  we  see  that 


e'^^^)^  (z  -  e  )    °  (5-7.15) 


We  now  examine  the  boundary  value  problem  for  the 
harmonic  function  Im{p(z){.   On  the  real  axis,  the  boundary 
values  are  0;  on  the  line  segment  of  slope  S   the  boundary 
values  are 

M(S^)    (■  _n    ) 

r\l\z   -   1^1    °  cos|M(S^)[7r  +  tan  ^S^]  j    (5-7-16) 

where  t^   Is  a  constant.   Thus,  by  employing  an  argument  similar 
to  that  used  above,  we  conclude  that  the  local  behavior  of 

A  . 

F(z},  as  far  as  any  singularity  Is  concerned.  Is 

M(S^) 
F(z)^  (z  -  ^J    °  (5-7.17) 

Thus,  In  the  neighborhood  of  the  "front"  or  "back"  of  the 
piston,  the  analytic  function 
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F(z)  =  P(z)e-^(^)  (5-7.18) 

will  be  regular. 

In  figure  21  we  have  sketched  the  local  geometry  of  the 
"top"  of  the  piston.   Here  we  have 


^o  =   ^^"    ^T=6S^)  (5-7.19) 


where  6  Is  the  (negative)  limiting  value  of  A  in  this  neigh- 
borhood, and  S  is  the  (negative)  limiting  slope  of  the  image 
of  the  piston.  Using  the  addition  formula  for  the  tangent, 
we  can  write 

A  -  tan"^6  +  tan"^S^  (5-7.20) 

To  discuss  the  local  behavior  of  the  harmonic  function 
Imlh(z)l  near  the  "top"  of  the  piston,  we  introduce  the 
harmonic  function 


N(S^)arg(z  -  ii]^)  (5-7.21) 


where 


tan"  5  +  tan  S 


N(S^)  =  -. (5-7.22) 

°      27r  +  tan  ^S^ 

This  function  is  chosen  to  have  the  same  local  boundary  values 
as  Im  lh(z)"l.   Thus,  in  the  neighborhood  of  the  "top"  of  the 
piston,  the  harmonic  function 
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Im[h(z)"^  -  N(S^)arg(z  -  1^^)         (5-7-23) 
will  assume  zero  boundary  values.   Thus,  locally. 


and 


h(z)  -  N(S^)log(z  -  iTi^)  (5-7.24) 


K/  ^  N(S^) 

.^^")'-  (z  -  mj    °  (5-7.25) 


To  proceed  to  discuss  the  boundary  value  problem  for 
F(z),  we  must  guarantee  that 


-1  <  N(Sq)  <  0  (5-7.26) 


Using  the  expansion  (4-3.15)5  and  recalling  that  in  the 
neighborhood  of  the  origin  the  mapping  from  the  {tt,   |-) -plane 
to  the  z-plane  Is  nearly  the  Identity  mapping,  we  can  write 

c^     2 
C  =  c[l  +  ^y^  n]  (5-7.27) 

Some  complicated  algebra  yields  the  approximation 

5  -  -  %  %[1  -  V  (4r  -  4r)  Ho]      (5-7.28) 
A  a    A 


so  that 

2 
A' 


tan  ^6  ft*  -  ^  n^  (5-7.29) 


Since  N(S^)  Is  Inherently  negative 
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-tan~^5  -  tan  ^S^  <  2Tr  +  tan'-^S^      (5-7.30) 


Using  (5-7-29)  we  can  write  (5-7-50)  in  the  form 

-tan'^S^  <  TT  -  ^  (^)  ri„  (5-7-31) 

A 

Since  r]     is  small ,  (5-7 .  31 )  is  always  satisfied,  as  is  (5-7.26) 

A 

Thus,  there  will  be  a  singularity  in  the  function  P(z). 

-h  f  z ) 
As  before,  multiplication  by  the  function  e   ^  '   will  remove 

this  singularity,  and  P(z)  will  therefore  be  regular.   Thus, 

the  solution  to  the  boundary  value  problem  will  always  be 

regular. 


5-8.   The  Normalization  of  the  Function  r 

Now  that  F(z)  has  been  constructed  we  can  use  equations 
(5_i|.19)  and  (5-4.20)  to  determine  g(C"^).   Since  F(z)  contains 
an  arbitrary  real  constant,  so  will  g(C  )•   If  we  use  the 
approximate  values 

&^  =  0  (5-8.1) 

5^  =  -1  (5-8.2) 

equations  (5-4.19)  and  (5-4.20)  can  be  written  as 
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^s  -    -^s^r  +  ^s^l  +  ^"  ^sS 


^s  =  ^s^  +  ^s^l  +  ^sS 


(5-8.3) 
(5-8.^) 


If  we  multiply  (5-8.3)  t)y  1^.  (5-8. '^)  by  \\    ,    and  subtract 


we  obtain 


's    's     r    ^  s 


+ 


nJC.g"      (5-8.5) 


's'  ^s 


where  we  have  used  the  fact  that 


^^s  ^<~-- 


(5-8.6) 


We  now  determine  the  value  of  the  arbitrary  real  constant 
that  appears  in  the  solution  by  satisfying  the  normalization 
condition  (5-3. '4-),  namely 


'O 


r(g,) 


dl  = 


IT 


(5-8.7) 


Equation  (5-3. 11)  serves  to  define  r(0  i^  terms  of  F(0  ^^'^ 


.+ 


iC)-      We  have 


r(0  = 


Ci 


pTD 


(5-8.8) 


The  variable  ^  is  real  and  here  lies  in  the  range 


c  <  C  <  C 


o 


(5-8.9) 


so  that  F(C)  in  the  denominator  of  (5-8.8)  is  identical  with 


F^(C) .      We    shall  write 


F^ia  =  R(o  +  r 


(5-8.10) 


where  J"*  is  the  undetermined  constant.   Equation  (5-8.5) 
gives 


g(^+)  = ^ — s 


so  that  we  can  write,  with  the  aid  of  (5-8.8) 


r(C) 


R(c)  +  r  +  il  -  % 

C3(R(0  +  r)(^3A-^  -  n^) 


(5-8.11) 


(5-8.12) 


Thus,  the  constant  "pis  to  be  determined  from  the  condition 


R(C)  +  P  +  ^s  -  % 


CC3(R(0  +  D^sA"-^  -  Tig) 


_ dC  =  I  (5-8.13) 


where 


s  =  s(0 


(5-8.14) 


Since  the  integration  in  (5-8.13)  extends  over  an 
interval  of  length  0(e  ),  and  the  right  hand  side  is  0(1), 
the  integrand  must  be  0(l/e  ).   This  means  that  F(z)  must 
have  nearly  a  zero  in  the  neighborhood  of  the  point  z  =  c. 
With  the  determination  of  the  value  of  the  constant  P  ,  all 
the  desired  properties  of  the  solution  are  fulfilled,  and 
all  the  functions  which  appear  in  the  solution  are  completely 
determined. 
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5-9-   Extension  to  Asymmetric  Pistons 

As  we  mentioned  earlier,  symmetry  about  the  (:g-)-axls 
Is  Inherent;  however,  we  can  allow  the  piston  cross-section 
to  be  asymmetric  about  the  (|-)-axls.   We  have  already 
Indicated  how  to  write  down  the  hyperbolic  wave  portion  of 
the  solution  for  such  a  piston. 

If  the  piston  Is  asymmetric  about  the  (|-)-axls  we  must 
Introduce  two  functions  r^^  and  r^.  ,  one  for  the  left  and  one 
for  the  right.   The  analytic  function  Is  obtained  by  solving 
a  half  plane  problem,  and  then  reflecting  across  the  (:r-)-axls. 
The  trouble  with  this  procedure  Is  that  we  are  no  longer  able 
to  make  F(z)  real  on  the  imaginary  axis.   Thus,  the  constant 
C   Introduced  In  (5-^.1)  will  not  necessarily  vanish. 

Since  we  are  matching  the  derivative  of  the  velocity 
rather  than  the  velocity  Itself,  the  vanishing  of  G   is 
essential  to  the  success  of  the  method.   Thus,  we  can  conclude 
that  there  will  be  solutions  for  an  asymmetric  piston  only 
If  the  geometry  is  such  that 


Im 


t 
f     (  2772 

^  "^2  ^  g ^  f(T)dT  =  0         (5-9.1) 

(t  -a  )(t  -c  ) 


This  condition  is  equivalent  to  a  condition  on  the  Green's 
function  of  the  domain  exterior  to  the  image  of  the  piston. 
It  is  reasonable  to  suspect  that  there  will  be  a  class  of 
asymmetric  pistons  for  which  conical  solutions  are  possible 
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6.   The  Solution  for  a  Symmetric  Piston  with  a  Partial 
Hyperbolic  Wave 

In  this  section  we  examine  the  solution  for  a  piston 
which  Is  rounded  at  the  "top"  and  "bottom."   In  this  case 
the  characteristics  drawn  from  the  cusped  portions  of  the 
slow  wave  envelope  (which  are  used  to  construct  the  hyper- 
bolic wave)  will  not  cover  the  whole  piston.   We  would 
expect  free  boundary  problems  In  the  neighborhood  of  the 
"top"  and  the  "bottom"  of  the  piston,  while  the  rest  of  the 
piston  (covered  by  the  hyperbolic  wave)  would  give  rise  to 
the  sort  of  boundary  value  problem  which  we  discussed  In  the 
last  section. 

In  this  section  we  consider  a  piston  cross-section  which 
Is  slightly  perturbed  from  the  circular  free  boundary  obtained 
earlier.   As  we  add  the  hyperbolic  wave,  we  seek  to  determine 
the  perturbations  to  the  free  boundary.   The  free  boundary  Is 
preserved,  to  first  order.  Indicating  that  mixed  solutions  of 
this  sort  are  possible.   It  should  be  noted  that  we  consider 
piston  expansion  velocities  so  small  that  the  shape  of  the 
piston  and  Its  Image  In  the  z-plane  are  coincident. 
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6-1.   Formulation  of  the  Problem 

In  this  subsection  we  denote  the  unperturbed  expansion 
speed  of  the  piston  by  r  ,  and  the  perturbation  to  that 
speed  by  er(0)  (see  figure  22).   We  write 

i{9)   =  (r^  +  £r(0))cos  9  (6-1.1) 

T^(e)  =  (r^  +  er(e))sin  6  (6-1.2) 

We  take  r(0)  to  be  unknown  in  the  angular  regions 


and 


)^  <  e  <  TT  -  9^  (6-1.3) 


TT  +  0^  <  0  <  27r  -  e^  (6-1.4) 


For  the  rest  of  the  piston  we  consider  r(0)  as  a  known 
function. 

In  section  3  we  obtained  the  expansion 

In  an  entirely  analogous  manner  we  obtain  the  expansion 

i=  (f)[l  +  i(|)+^(|)'  +  0(r5)j        (6-1.6) 

x    y 
Since  the  point  {-^,   -^)    shown  in  figure  22  is  the  point  of 

intersection  of  the  last  real  characteristic  from  the  slow 
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wave  envelope,  and  the  piston,  we  can  match  the  slope  of 
the  piston  at  that  point  with  the  slope  of  the  hyperbolic 
wave.   The  slope  of  the  hyperbolic  wave  is 

"(f-f)=^=-|(f)t+|(|)  +  °(=^o)}   (6-1.7) 

Since  the  piston  is  nearly  circular,  its  slope  is  approxi- 
mately given  by 

S(f,  f)  =  -(f)/{|)  (6-1.8) 


We  now  set 


^O   ^On      ,^0   ^O. 


The  approximate  solution  of  (6-1.9)  is 

2 


X            r 

(^)«r„(l- 

2 
r 

o 

(6-1.10) 


(6-1.11) 


IT 


Thus,  we  see  that  0   is  very  close  to  ^-^   On  the  small 
intervals  at  the  "top"  and  "bottom"  of  the  piston  we  solve  a 
free  boundary  problem.   On  the  rest  of  the  piston  we  solve  a 
boundary  value  problem  for  both  the  elliptic  solution  and  the 
hyperbolic  wave.   A  general  existence  theorem  is  very  diffi- 
cult to  provide;  the  analytic  function  f(z)  which  results 
from  the  various  portions  of  the  solution  must  be  shown  to  be 
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one  analytic  function  In  the  whole  z-plane,  excluding  the 
Image  of  the  piston. 


6-2,   The  Solution  of  the  Mixed  Problem 

Using  the  notation  of  section  5>  we  write  the  function 

F(z)  as 

2 

F(z)  =  ^  +  ef(z)  (6-2.1) 

z 

Using  (6-1.1)  and  (6-1.2)  for  the  regions  where  r(0)  Is  known 
we  find  that  the  analogue  of  equation  (5— 4-.23)  Is  the  linear 
relationship 

(^0+A^g)F^(^,Ti)  +  (e^-Ari0)P^(^,r,)  =  Hg  -A^g     (6-2.2) 

The  order  one  terms  of  this  equation  constitute  an  Identity; 
the  order  e  terms  yield 

a{e)f^{9)   -  p(e)f.(0)  =  y(9)  (6-2.5) 

where 

a(0)  =  cos  9  -  A  sin  9  (6-2.4) 

P(0)  =  sin  9  +   A   cos  9  (6-2.5) 

t(9)  =  ■§-   [r(0)(cos  e+A  sin  9)   +rn(e)(sln  0  -A  cos  9)] 
o 

(6-2.6) 
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Prom  section  3  we  know  that  in  the  two  Intervals 
0     <   9  <   IT  -   9     and  t  +   9     <0<2ir-0   (where  the  radial 
perturbation  r(0)  Is  unknown)  we  must  satisfy  the  two  relations 

fI  +  fI  =   1  (6-2.7) 

where  S  is  the  slope  of  the  image  of  the  piston  in  the  z-plane. 
Together  these  two  equations  imply  that 


nQ{0)F^{9)    =  ^g(0)  +  iQ{9)F^{9)  (6-2.9) 

The  order  one  terms  again  constitute  an  identity;  the  first 
order  terms  in  the  parameter  e  yield 

a(0)f^(0)  -^(0)f.(0)  =  7(0)  (6-2.10) 
where 

a(0)  -  -  sin  0  (6-2.11) 

p(0)  =  cos  0  (6-2.12) 

^(9)  =  ^  (r(0)sin  0  -  rg(0)cos  0)  (6-2.13) 

o 

In  the  two  small  angular  regions  where  these  relations 
are  valid  we  can  set 

cos  0  =  0  (6-2.14) 

sin  0  =  ±  1  (6-2.15) 
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where  the  plus  sign  in  (6-2.15)  refers  to  the  "top"  of  the 
piston,  and  the  minus  sign  refers  to  the  "bottom."   Thus, 
at  the  "top"  of  the  piston  (6-2.10)  gives 

?r(0)  =  ^4^  (6-2.16) 

while  at  the  "bottom"  of  the  piston  this  equation  gives 

J  (e)  =  -  2ZM  (g_2.i^) 

o 
Since  we  want  F(z)  to  be  real  on  the  real  axis,  we  must 

A 

have  that  f   is  an  even  function,  upon  reflection  in  the 
r 

real  axis.   Thus,  we  conclude  that 

r(0)  =  0  (6-2.18) 


m 


(0Q  <  9  <  TT  -  9^)  U  (tt  +  9^  <  9  <  27r  -  9^)  (6-2.19) 


Thus,  we  now  know  the  entire  cross-section  of  the  piston, 
and  can  employ  the  methods  of  the  last  section  to  complete 
the  solution.   This  result  indicates  that  there  are  solutions 
to  mixed  problems  of  this  sort. 
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7-   Solution  for  a  Thin  Piston  with  Finite  Expansion  Velocity 

In  this  section  we  obtain  the  solution  for  a  symmetric 
piston  whose  surface  Is  nearly  parallel  to  the  (tt) -axis .      As 
would  be  expected.  In  light  of  our  earlier  remarks,  the 
linearization  Is  valid  In  this  case.  Inasmuch  as  the  pertur- 
bations to  the  flow  variables  are  all  0(e),  where  here  e  Is 
the  maximum  of  the  (^) -component  of  the  piston  velocity.   We 
take  the  maximum  (:p) -component  of  the  velocity  to  be  0(l)  but 
less  than  c  In  magnitude  (see  figure  23).   Finally,  we 
assume  that  the  entire  surface  of  the  piston  Is  covered  by 
the  hyperbolic  wave. 

Since  the  neighborhood  of  the  origin,  and  the  (^)-axls 
are  preserved  under  the  mapping  from  the  (ir,  p) -plane  to 
the  z-plane,  we  can  take  the  Image  of  the  piston  to  be  Iden- 
tical with  the  piston  Itself. 

In  the  first  subsection  we  discuss  the  boundary  conditions 
on  the  surface  of  the  piston.   The  conditions  we  obtain  are 
slightly  different  from  those  obtained  In  section  one.   In 
particular,  we  take  B-n  =  0  on  the  surface  of  the  piston, 
so  that  this  surface  Is  flux  preserving. 

The  second  subsection  discusses  the  reduction  of  the 
problem  to  the  sort  discussed  In  section  five.   In  the  last 
subsection,  we  treat  the  piston  Image  as  a  silt  In  the 
z-plane,  and  write  down  an  explicit  formula  for  the  analytic 
function  F(z).   We  do  this  by  employing  the  Plemelj  Integral 
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f  orm-ulas . 

It  Is  worth  noting  that  there  appears  to  be  no  solution 
for  a  thin  piston  with  expansion  speed  greater  than  c.   In 
this  case  the  hyperbolic  wave  will  not  cover  the  whole  piston, 
and  the  remaining  portion  of  the  piston  does  not  resemble 
the  circular  free  boundary  obtained  In  section  3- 


7-1 .   Boundary  Conditions  on  the  Piston 

On  the  surface  of  the  piston  we  require  that 

u-n  =  U  -n  (7-1-1) 

B-n  =  B  -n  (7-1.2) 

p  \  '  / 

and 

B-n  =  0  (7-1-5) 

This  last  condition  Insures  that  the  magnetic  field  Is 
tangential  to  the  surface  of  the  piston,  and  that  the  surface 
Is  flux  preserving.   With  B  =  0  It  Is  no  longer  appropriate 
to  prescribe  [E^]  =0,  [6],  since  the  former  condition  corre- 
sponds to  the  loss  of  a  characteristic  cone. 

The  fields  In  the  piston  are  determined  by  solving  the 
pre-Maxwell  equations  (without  displacement  current)  with 
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the  boundary  conditions 


B  -n  =  0  (7-1.^) 

P 


E  +  U  X  B  =  0  (7-1-5) 

P    P    P 

This  is  a  well-posed  problem  for  the  pre-Maxwell  equations  [3] 

If  we  denote  the  normal  vector  to  the  surface  of  the 
piston  by 

n  =  (n  ,n  )  (7-1-6) 

^  X   y^  \ I     / 

we  see  from  figure  23  that 

n^  =  0(e)  (7-1.7) 

Hy  =   0(1)  (7-1-8) 

Using  the  linearization  (1-1.8)  we  can  write  (7-1-1)  and 
(7-1.2)  in  the  form 


un  +  vn  =  U,n  +  U^n  (7-1-9) 

X      y     1  X     2  y  \  i    ^  / 


Bn  +Bn  =-n  (7-1.10) 

X  X    y  y      X  \  r      / 

Since 

—  =  -  S  (7-1-11 

n 

y 

where  S  is  the  slope  of  the  cross-section  of  the  piston  we 


-  109  - 


finally  have  the  boundary  conditions 


V  =  -  SU-j_  +  Ug 


B  =  +  S 

y 


(7-1.12) 
(7-1.13) 


7-2 .   Solution  of  the  Boundary  "Value  Problem 


We  adopt  the  formalism  and  notation  of  section  five, 
and  write 


C, 


f(2)  = 


(      ro 

2  I    ^r(^) 

TT 


1-^^  dl  [  F(z) 


(7-2.1) 


In  this  case,  however,  we  normalize  the  function  r  by 
requiring 


'O 


er(^)d^  =.  I 


(7-2.2) 


and  on  the  surface  of  the  piston  set 


f(z)  ^^^ 


2         2 
c  -  z 


(7-2.3) 


We  then  satisfy  the  derivatives  of  (7-1.12)  and  (7-1.13), 
namely. 
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dv 
ds 


^  d?  ^C^ 


(7-2. t) 


ds 


_d_ 

ds  ^^ 


(r^) 


(7-2.5) 


Straightforward  algebra  yields  the  two  equations 


dv 
ds 


^  7~2   2  £2v2    ,,.2  2   ^  ^r 
(c  +T1  -I  )   +  ^^  11 

(cW-^^)^3  -  2|titi^ 

"^  ^  (    ^,  2  ,2^2  ^  .,2^2   ^  ^1  +  ^sS 
(c  +11  -e,  )  +  4^  n 


(7-2.6) 


dB^^     -2^^nl  +Ti(cW-^^)?.  -  ^(c^+t1^-|^)ti  -2^Ti2ri 

y  _  r 2 2_ £ £.  1 T? 

:^   2vr  /  2   2  -.2x2      ^""5  -l-^- 


ds 


[■ 


(r+n")[(c"+Ti"-r)"  +  ^Frn 


-2^n^lg  -  ^(c^+n^-4^)^g  -  il(c^+il  -l^)Tlg  +2^^riTi^ 

+  [  TT ^ ^ ?; ?5— ^ T^—TT ]  F. 

(r+ri2)[(cW-r)^  +  ^erf]  ^ 


-  ^s^f) 


(7-2.7) 


An  examination  of  the  relative  parity  of  these  two  equations 
shows  that  the  relative  parity  of  the  real  and  imaginary 
parts  of  the  analytic  function  P(z)  are 


+ 


P, 


+ 


+ 


+ 


P.: 


+ 


+ 


while  the  relative  parity  of  the  function  g  is 
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+ 
+ 


This  is  just  as  we  found  In  section  five.   Thus  we  can  again 
consider  a  quarter  plane  boundary  value  problem,  and  extend 
the  resulting  function  to  the  whole  plane  with  two  applica- 
tions of  the  reflection  principle. 

If  we  now  drop  lower  order  terms  In  (7-2.6)  and  (7-2.7) 
we  can  write  (7-2.^)  and  (7-2.5)  In  the  form 

-  -J-  (7-2.9) 

We  now  eliminate  g  from  these  two  equations  to  obtain  a 
linear  equation  connecting  the  boundary  values  of  the  real 
and  Imaginary  parts  of  F.   We  write  this  relation  as 

aF^  +  PF^  =  Y  (7-2.10) 

where 

_  [-2^Ti^^-(c2-^2)^^](^2^^2)^^^2^2^^^_^^^^_^^^^^^2_^2^^ 
a ^ — ^_^ 


(r+T  )(c  -e")' 

(7-2.11) 
a,  § c^      o        n      o        ^ ^ (7-2.12) 

(r+n^)(c2-^2) 


-  112  - 


(^  +ri  )(c  -^  ) 
and 

Y  =  (^-C)  ^  (^)  (7-2.14) 

s 


We  take  the  piston  to  be  smooth,  so  that 

dl  (^)  -  0(e)  (7-2.15) 

s 

Then  we  can  calculate  that 

a  =  0(e)  (7-2.16) 

P  =  0(1)  (7-2.17) 

T  -  0(e)  (7-2.18) 

One  might  hope  that  further  simplification  could  be  obtained 
by  setting  a  =  0.   However,  this  reduction  sufficiently 
changes  the  character  of  the  problem  so  that  no  solution  is 
possible . 

We  have  discussed  boundary  value  problems  of  the  sort 
(7-2.10)  in  section  five.   To  be  able  to  carry  out  the 
solution,  we  must  insure  that 


,    -1  /Ctfs 
tan 


■f3(sj)  (7-2-19) 

remains  on  one  branch  of  the  arctangent  function.   We  make 
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the  approximation  C  =  c,  and  find  from  (7-2.12)  and 
(7-2.13)  that 

^    2|2ri-(c+e)[^Vc(n-|S)] 


(7-2.20) 


P        (c+^)[^(c-n-il^] 
The  denominator  of  this  expression  has  a  zero  at 

n  =  +  y^{c-i)  (7-2.21) 

The  numerator  can  be  written  in  the  form 

(c-^)[e(c+e)S  -  Tl(c  +  2^)]  (7-2.22) 

This  expression  has  a  zero  only  for  S  >  0.   In  the  first 
quadrant  S  <  0,  so  that  the  numerator  of  (7-2.20)  has  a 
zero  only  when  t]  =  0  and  then  only  if  S  =  0.   Thus,  if  we 
take  the  branch 


0  <  tan  ^0  <  TT  (7-2.23) 


the  boundary  values  we  specify  for  the  imaginary  part  of 
the  auxiliary  analytic  function  h(z)  all  lie  on  one  branch 
of  the  arctangent.   Thus,  the  problem  reduces  to  the  problem 
considered  in  section  five,  except  here,  on  the  piston 

|f(z)|  =   0(e)  (7-2.24) 

since  the  final  boundary  data  are  all  0(e).   Thus  in  this 
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case  the  linearized  problem  appears  to  be  a  faithful  model 
of  the  nonlinear  problem. 

Our  uniqueness  theorem  no  longer  holds,  since  In  this 
case 


u  /=  U-^  (7-2.25) 


on  the  piston.   The  right  hand  side  of  equation  (1-6.6)  is 

2 
no  longer  zero;  it  is  a  quantity  of  0(£  ). 


7-3-   An  Explicit  Determination  of  the  Function  F 

Referring  to  the  method  of  section  five  to  solve  the 
sort  of  boundary  value  problem  posed  by  (7-2.10),  we  can 
write  down  an  explicit  formula  for  the  analytic  function 
F(z).   We  treat  the  piston  image  as  a  slit  in  the  z -plane, 
and  employ  the  Plemelj  formulas. 

We  write 


and  note  that 


A(U  =  tan-^(|^)  (7-3.1) 


A(-e)  =  -  A(^)  (7-3.2) 
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We  denote  the  maximum  of  the  (:^) -component  of  the  piston 
velocity  by  [x,    and  then  write  the  analytic  function  h(z)  as 


h(z)  =i 


IT 


+\^ 


-^x 


T  -  Z 


dT 


(7-5.3) 


On  the  real  axis,  In  the  Interval  -ix  <  ^  <  +ix  we  have 

P  I   A(t' 


h  (4,0)  = 


dT 


(7-3.^) 


The  boundary  values  of  Im  |  P(z)|  on  the  silt  from  -a  to  a 

h^(4.0) 


are 


k(|)  =  sgn(4)Y(4)e 


(7-3.5) 


We  can  then  write 


F(z)  = 


T 


IJ- 


-l-L 


^1:^1  dT  +  r 


T  -  Z 


(7-3.6) 


where  p  Is  an  arbitrary  real  constant.   The  function  F(z) 
Is  then  given  by 


F(z)  =  F(z)e 


-h(z) 


(7-3.7) 


It  Is  easily  verified  that  this  function  satisfies  all  the 
desired  conditions.   The  constant  P*  Is  determined  from  the 
normalization  condition  (7-2.2). 
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Finally,  it  is  an  easy  matter  to  estimate  the  behavior 

of  F(z)  near  the  points  ±  |i.  and  0,  where  it  may  possess 

singularities.   This  is  a  question  of  the  behavior  of  Cauchy 

type  Integrals  near  the  endpoints  of  the  Interval  of 

Integration.   (To  discuss  the  point  0  we  write 

+\i        o    +\i 
J     =   J     +   J        ,  so  that  zero  becomes  an  endpoint . )   To 
-\x        -\i        o 
discuss  the  behavior  of 


t  -z 


dt 


(7-3.8) 


in  the  neighborhood  of  the  point   a,  when  i/'(a)  <  oo^  ^e  write 


C(.).^    iltl^llai  ,,  ,  |M 


dt 
t-z 


(7-3.9) 


or 


(^)  =  ^o(-)  +  1^  l°«(f^' 


(7-3.10) 


where 


b 


2it1 


•^(t)  --^(a)  ^^ 
t-z 


a 


(7-3.11) 


will  be  analytic  in  the  neighborhood  of  z  =  a.   We  take  a 
definite  branch  of  the  logarithm,  which  is  single-valued 
in  the  plane  slit  from  a  to  b.   Then  we  can  write 


C(z)  .  c^(.)  .  |M  log(^) 


(7-3.12) 
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where  C  (z)  Is  analytic  In  the  neighborhood  of  z  =  a.   A 
similar  arg-ument  yields  the  estimate 

C{z)  =  C,(z)  -  |M  log(^)  (7-3.13) 


where  C,  (z)  is  analytic  in  the  neighborhood  of  z  =  b , 
In  the  case  that 


^(t)  =  —^ (7-3.1^) 


(t  -a) 


a 


where  T^^ft)  is  regular  at  z  =  a,  and  0  <  a  <  1,  an  examination 

n 


of  the  function 


n(z)  = 


b 

^1 '  (7-3.15) 

(t-a)^(t-z) 


shows  that  C(z)  has  a  singularity  of  the  form 

(7-3.16) 


(z-a) 


a 


near  the  point  a . 

We  will  now  examine  the  function  h(z)  near  z  =  0  and 
z   -  +\i    .      If  the  piston  image  intersects  the  positive  ^-axis 
with  a  slope  S  ,  the  limiting  value  of  'K{i) ,    from  equations 
(7-3.1)  and  (7-2.20)  is 

A^  =  tan-^S^  (7-3.17) 


-  Hi 


Thus,  near  z  =  (x,  h(z)  has  the  local  behavior 


h(ii)  ~  log[  (z-[i)   °  ]  (7-3.18) 


where 


tan  S 
M(S^)  ^   -   ( ^)  (7-3.19) 


By  our  choice  of  branch  for  the  arctangent,  0  <  M(S  )  <  1 


Thus,  from  (7-3.5)  we  have 


k(h)  ^  (n-c)  ^  (S^)  (7-3.20) 


so  that 


where 


F(^l)  --  N(S^)log(z-^L)  +  P  (7-3.21) 


N(S^)  =  -  4^  (7-3.22) 


IT 


Thus,  the  function 


P(z)  =  F(z)e"^^^)  (7-3.23) 


has  a  singularity  of  the  form 


-M(S  ) 
(z-ii)  °  (7-3.2^) 


which  will  always  be  Integrable,  and  thus  cause  no  singularity 
In  the  flow. 
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At  the  point  z  =  V,  where  the  piston  Image  Intersects 
the  positive  n-axls,  A  tends  to  the  limiting  value 


-1/C 


A^  =  tan--^(^)  (7-5.25) 


Near  z  =  v,    h(z)  has  the  local  behavior 

h(v)  ^   log[(z-v)^]  (7-3.26) 

where  here 

M  =  -  tan-^(c/v)  (7-3.27) 

IT 

Here  k(v)  Is  not  finite,  so  we  must  use  our  second  result 
(7-3-16)  which  gives 

F(v)  ^  (z-v)^  (7-3.28) 

Thus,  the  function 

F(z)  =  P(z)e"^^^^  (7-3.29) 

will  be  regular  at  z  =  v. 

Summing  up,  the  function  F  will  have  weak  Integrable 
singularities  at  the  "tips"  of  the  piston,  but  the  resulting 
flow  will  be  regular. 
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Figure  Captions 

1.  The  normal  speeds  locus 

2.  Propagation  of  plane  waves 

3.  The  characteristic  conoid 

4.  A  t  =  constant  section  of  the  characteristic  conoid 
5'  Tangents  to  the  characteristic  conoid 

6.  Tangents  to  the  characteristic  conoid 

7.  Placement  of  the  piston  cross-section 

8.  The  z -plane 

9.  The  image  of  the  piston  in  the  z -plane 

10.  Definition  of  arc  length  on  the  Image  of  the  piston 

11.  Mapping  of  the  z-plane  into  the  ^ -plane 

12.  Geometry  of  the  mapping 

15.  Mapping  of  the  z-plane  into  the  w-plane 

ih .  The  real  characteristics 

15.  Geometry  of  the  hyperbolic  wave 

16.  Details  of  the  geometry  of  the  hyperbolic  wave 

17.  The  z-plane  and  the  image  of  the  piston 

18.  Quarter-plane  problem  for  Im|h(z)] 

19.  Quarter-plane  problem  for  Im^F(z)'l 

20.  Local  boundary  value  problem  for  Im^h(z)] 

21.  Local  boundary  value  problem  for  Im^h(z)\ 

22.  Piston  for  the  mixed  problem 

23.  Finite-size  piston 
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Figure  3 
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Figure  8 
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Figure   12 


Figure  13 


Figure    14 
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Figure    15 
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Figure  17 
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Figure  21 
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